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Abstract 

Within  the  last  decade,  many  new  technologies  have  significantly  changed  the  face 
of  private  astronomy.  Developments  such  as  inexpensive  but  high-quality  sensors,  rapid 
personal  computing,  and  easy  networking  inspire  a  reexamination  of  an  old  problem: 
how  practical  is  it  to  develop  initial  orbit  estimates  for  Low  Earth  Orbiting  (LEO)  satel¬ 
lites  using  optical  tracking?  This  paper  documents  the  design  and  implementation  of  a 
commercial  telescope  system  used  to  answer  precisely  that  question.  This  analysis  deter¬ 
mined  there  are  some  challenging  barriers  to  successful  single-site  orbit  determination, 
but  it  is  possible  given  the  right  conditions.  Considering  the  low  cost  and  small  sup¬ 
port  footprint  of  such  systems,  they  could  provide  excellent  support  to  Space  Situational 
Awareness  (SSA)  missions  or  satellite  tracking  operations  in  general. 
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Initial  Determination  of  Low  Earth  Orbits  Using 


Commercial  Telescopes 


I.  Problem  Statement 


RAPID  advances  in  the  quality  of  electronics,  combined  with  equally  dramatic  im¬ 
provements  in  cost  and  availability,  are  revolutionizing  private  astronomy.  Unpar¬ 
alleled  access  to  quality  equipment,  rapid  personal  computing,  and  extensive  community 
support  enable  nearly  anyone  to  achieve  feats  in  their  backyard  that  required  an  obser¬ 
vatory  twenty  years  ago.  Semi-professional  astronomers  and  programmers  continually 
develop  novel,  inexpensive  methods  to  defeat  complex  engineering  challenges. 

One  such  challenge  is  optically  tracking  satellites  to  determine  their  orbits.  There 
are  long-standing  solutions,  but  this  project  approaches  it  with  refreshed  interest.  Pri¬ 
mary  motivations  include: 

•  Space  Situational  Awareness  (SSA):  Commercial  systems  are  inexpensive,  mobile, 
and  easily  supported:  all  factors  that  compensate  for  limitations  in  capability. 
There  are  always  new  opportunities  to  use  them  for  surveillance  and  debris  moni¬ 
toring. 

•  State-of-the  Art  Survey:  This  study  offers  a  baseline  explanation  of  methods  used 
by  semi-professional  satellite  observers  today.  Hopefully,  it  will  serve  as  a  reference 
for  future  researchers  and  motivate  them  to  pursue  additional  work  in  this  held. 

•  Research  Testbed:  By  its  conclusion,  this  project  integrated  the  hardware  and  soft¬ 
ware  required  to  operate  a  basic  optical  satellite  tracking  program.  Now,  students 
may  use  it  to  support  work  in  sensors,  image  processing,  orbit  determination,  and 
many  other  holds.  It  also  allows  AFIT  students  to  gain  hands-on  experience  with 
classroom  concepts. 
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Tracking  artificial  satellites  is  a  pastime  as  old  as  the  Space  Age  itself.  Like  many 
other  subdisciplines  of  astronomy,  this  held  benefits  greatly  from  recent  advances.  This 
project  examines  how  modern  equipment  is  used  to  track  Low  Earth  Orbit  (LEO)  satel¬ 
lites  in  order  to  determine  their  orbits. 

Successful  determination  of  any  body’s  orbit  requires  accurate  measurements  in 
both  space  and  time.  For  millennia,  astronomers  had  no  way  of  estimating  how  far  away 
stars  and  planets  were,  so  the  methods  they  developed  for  predicting  their  positions 
relied  only  on  time  and  relative  angular  measurements.  Great  minds  of  the  day  (namely 
Gauss  and  Laplace)  developed  very  robust  routines  for  calculating  orbits  with  such  data. 
Today,  these  methods  are  collectively  called  “angles  only”  techniques. 

When  artificial  earth  satellites  were  first  launched  in  the  late  1950s,  they  brought 
with  them  a  pressing  need  for  accurate  orbital  measurements.  Without  them,  it  would 
be  impossible  to  keep  track  of  a  launched  satellite.  Astronomers  implemented  familiar 
angles  only  methods  to  monitor  these  new  celestial  bodies,  yet  they  would  face  unique 
new  challenges.  Chapter  II  describes  the  origins  of  optical  satellite  tracking  and  enu¬ 
merates  key  data  needs  for  successful  orbit  determination.  It  also  describes  how,  using 
fundamental  principles  established  decades  ago,  it  is  possible  to  apply  modern  comput¬ 
ing,  imaging,  precision  navigation,  and  timing  technology  to  produce  effective  results 
at  low  cost.  Chapter  II  concludes  with  a  brief  overview  of  the  hardware  and  software 
selected  for  this  project. 

Chapter  III  discusses  the  prerequisite  task  of  identifying  opportunities  for  visual 
satellite  tracking.  First,  the  North  American  Aerospace  Defense  Command  (NORAD) 
satellite  catalog  is  discussed,  followed  by  an  overview  of  the  SGP4  algorithm  used  to 
extract  orbits  contained  in  the  catalog.  Then,  key  transformations  between  local  and 
inertial  frames  are  described  in  detail.  Finally,  satellite  brightness  models  are  applied. 
Once  these  steps  are  complete,  it  is  possible  to  pursue  observations  with  a  high  certainty 
of  success. 

Chapter  IV  explains  how  collected  data  is  processed  to  produce  local  angular  mea¬ 
surements  of  satellite  overflights.  Then,  these  measurements  are  transformed  to  an  iner¬ 
tial  reference  frame  in  preparation  for  initial  orbit  determination. 
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The  results  of  system  calibration  and  early  observations  are  presented  in  Chapter 
V.  A  case  study  is  presented  that  reveals  this  system  is  capable  of  producing  useful  initial 
orbit  estimates,  given  the  right  conditions.  There  are  numerous  theoretical  and  practical 
concerns  that  complicate  tracking  with  a  single  telescope  from  a  single  site.  Considering 
the  low  cost  of  systems  like  this,  however,  these  complications  may  become  irrelevant  if 
multiple  sensors  are  employed. 

Chapter  VI  summarizes  the  project  and  expands  on  the  initial  conclusions  deter¬ 
mined  in  this  course  of  study.  It  identifies  specific  research  areas  that  deserve  further 
analysis,  in  order  to  both  deepen  understanding  of  optical  tracking  and  build  a  founda¬ 
tion  future  tracking  networks  can  rest  on. 
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II.  Background 


From  the  ground,  an  observer  can  spot  sunlight  reflecting  off  of  a  satellite  under  the 
right  conditions,  much  like  sunlight  reflects  off  the  moon.  This  means  a  satellite 
can  appear  as  bright  as  a  star,  with  the  notable  exception  that  it  moves  much  faster 
across  the  night  sky. 


In  October  of  1957,  this  effect  caused  great  concern  in  the  United  States  as  the 
world’s  first  artificial  satellite,  Sputnik  1,  could  easily  be  seen  flying  methodically  over¬ 
head.  The  shiny  metal  sphere’s  polish  went  beyond  propaganda  -  its  mirror-like  surface 
aided  telescope  tracking  [ Smithsonian ,  2008].  The  probe’s  radio  beacon  had  a  limited 
lifespan,  and  only  one  radar  in  England  was  capable  of  tracking  the  relatively  large  rocket 
body  that  remained  in  space,  not  Sputnik  itself  [BBC,  2007]. 


Figure  2.1:  Sputnik  1  [NASA,  2007] 


Shortly  before  Sputnik’s  launch,  the  Smithsonian  Institution  organized  Operation 
Moonwatch.  Its  international  observer  corps  first  tracked  Sputnik,  then  other  satel¬ 
lites  over  many  years  that  followed.  Volunteers  used  arrays  of  very  simple  instruments  to 
record  the  time  and  place  a  target  satellite  passed  a  given  observing  site.  These  measure¬ 
ments  were  used  to  calculate  satellite  orbits  and  also  determine  geophysical  properties 
of  the  Earth  and  its  atmosphere.1 

1For  a  thorough  and  delightfully  campy  history  of  the  early  days  of  satellite  tracking,  refer  to  [Engle 
and  Drummond,  1965]. 
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AN  AUSTRALIAN  Moonwatch  Team  is  part  of  a  world  wide  group  of  volunteer  observers 
who  have  made  over  100,000  valuable  observations  of  satellites.  Moonwatch  observers 
can  compete  with  photographic  and  radar  stations  especially  on  very  low  altitude  fast 
moving  objects.  They  have  proven  to  be  a  necessary  part  of  the  satellite  tracking  program. 


Figure  2.2:  Australian  Project  Moonwatch  Volunteers  ( Photo  Courtesy  the  Harvard- 
Smithsonian  Center  for  Astrophysics) 

As  the  U.S.  government  developed  a  comprehensive  satellite  tracking  network,  how¬ 
ever,  the  role  of  optical  tracking  changed  dramatically.  Large  sky-scanning  radars  were 
developed  that  could  find  and  measure  objects  in  LEO  regardless  of  sky  conditions.  Ad¬ 
ditionally,  radars  can  determine  range  and  instantaneous  change  in  range  (range-rate) 
to  the  target.  Because  of  this,  telescopes’  missions  narrowed.  Today,  the  Ground-Based 
Electro-Optical  Deep  Space  Surveillance  (GEODSS)  system  is  used  to  observe  objects 
beyond  radar  range,  developing  orbits  for  objects  10,000  to  45,000  kilometers  from  Earth. 
GEODSS  telescopes  have  a  one  meter  diameter  and  use  a  Charge-Coupled  Device  (CCD) 
camera  to  detect  objects  10,000  times  dimmer  than  visible  with  the  naked  eye  [USAF, 
2006].  Undoubtedly,  GEODSS  is  many  times  more  capable  than  its  predecessors,  and 
represents  the  technological  peak  of  optical  satellite  tracking. 

Two  other  government-sponsored  systems  are  noteworthy  as  well,  because  they  use 
commercial  equipment  to  track  and  analyze  satellites.  The  United  States  uses  the  Raven 
system  to  support  its  space  tracking  efforts,  including  satellite  and  debris  characterization 
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and  tracking  in  all  orbital  regimes  [. Kervin  et  al.\.  The  Canadian  Satellite  Tracking  and 
Orbit  Research  (CASTOR)  has  stated  objectives  to  bring  amateurs  and  professionals 
together,  and  performs  varied  missions  as  well.  In  honor  of  Sputnik’s  50th  anniversary, 
CASTOR  tracked  over  2000  distinct  objects  in  2007  [Earl,  2008].  Both  systems  are  now 
in  service  for  over  a  decade. 

Whether  using  a  simple  or  sophisticated  telescope,  their  fundamental  purpose  in 
orbit  determination  is  the  same:  accurately  measure  the  apparent  position  of  a  target  at 
a  specific  time.  The  remainder  of  this  chapter  discusses  bedrock  concepts  that  help  reach 
that  seemingly  simple  goal.  Since  this  project  deals  exclusively  with  observing  satellites 
in  the  visible  range,  Section  2.1  describes  how  astronomers  describe  the  brightness  of  a 
celestial  body.  Section  2.2  describes  what  theoretical  parameters  are  required  to  develop 
orbits  using  telescope  data,  followed  by  a  section  on  major  sources  of  potential  error. 
Finally,  Section  2.3  provides  a  brief  overview  of  modern  techniques  used  to  eliminate 
major  observational  errors  and  explains  the  equipment  used  in  this  study. 

2.1  The  Visual  Magnitude  Scale 

Quantifying  the  brightness  of  nighttime  objects  is  hardly  a  recent  pursuit.  The 
ancient  Greek  astronomer  Hipparchus  developed  a  catalog  of  stars’  intensities  by  subjec¬ 
tively  placing  any  one  he  could  see  into  one  of  six  categories.  Ptolemy  would  continue 
the  tradition  of  referring  to  the  brightest  stars  as  first  magnitude,  whereas  sixth  mag¬ 
nitude  stars  were  barely  perceptible  [ Kennon ,  1948].  Despite  two  millennia  of  scientific 
advance,  this  terminology  and  the  concept  of  apparent  magnitude  remains  embedded  in 
the  lingua  franca  of  astronomy.  The  following  concepts  are  critical: 

•  Higher  is  Lower:  Unlike  most  scientific  scales,  brighter  objects  have  numerically 
lower  apparent  magnitudes,  extending  into  the  negative  range. 

•  Logarithmic  Scale:  A  decrease  of  one  in  apparent  magnitude  (say,  5  to  4)  corre¬ 
sponds  to  a  v^lOO  ~  2.512  multiplication  in  brightness.  This  ratio  was  proposed 
by  Norman  Pogson  as  a  standard  in  1856  [ Pogson ,  1856].  To  express  any  scalar 
multiple  in  apparent  magnitude  (Mx)  as  a  common  logarithm,  the  following  change 
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of  base  is  used: 


log(Wa>)(M,)  =  i^a^b  =  2.51og10(M„)  (2.1) 

The  apparent  magnitude’s  logarithmic  scale  is  well  suited  for  the  human  eye’s 
logarithmic  visual  response  curve.  Other  sensors,  such  as  him  or  CCD  cameras, 
respond  differently.  Refer  to  [Rees,  2001]  for  an  introduction  to  remote  sensing 
methods  that  take  sensor  performance  into  account. 

•  Absolute  Magnitude:  A  reference  brightness  called  absolute  magnitude  defines  an 
object’s  brightness  if  observed  at  a  standard  distance  and  orientation.  Depending 
on  the  desired  correction,  an  apparent  magnitude  may  be  predicted  for  a  given 
geometry [Pogson,  1856;  Meeus ,  1998]. 

Pogson  is  credited  with  creating  the  modern  visual  magnitude  scale.  During  his 
time,  astronomers  were  searching  for  a  suitable  model  to  predict  the  brightness  of  aster¬ 
oids.  Pogson  claimed  his  formula  could  accurately  match  observed  trends,  and  that  any 
errors  would  be  constant.  He  predicted  this  constant  error  would  be  caused  only  by  his 
mis-estimation  of  each  body’s  absolute  magnitude;  he  then  boldly  notes  “but  this  I  do 
not  anticipate”  [ Pogson ,  1856].  This  philosophy  is  simple:  pick  a  model  that  matches 
trends  and  adjust  the  bias  (in  this  case  absolute  magnitude)  to  match  observations. 

A  century  later,  F.L.  Whipple  and  J.A.  Hynek  calculated  apparent  magnitude 
estimates  for  orbiting  satellites  as  they  designed  the  United  States’  first  tracking  telescope 
network.  They  presented  their  working  assumptions  to  the  Institute  of  Radio  Engineers: 

Calculations  show  that  at  a  zenithal  distance  of  about  200  miles  in  twilight  a 
20-inch  sphere  with  albedo  0.6  would  have  a  photographic  magnitude  of  6.3 
and  a  visual  [apparent]  magnitude  of  5.7.  [Whipple  and  Hynek,  1956] 

The  exact  calculation  and  assumptions  they  used  are  undocumented,  but  Chapter  III, 
Section  3.3  will  describe  how  methods  much  like  those  Pogson  and  Whipple  employed 
are  still  in  use  today. 
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2.2  Angles  Only  Orbit  Determination 

Orbit  solutions  that  do  not  use  target  range  data  are  collectively  called  angles  only 
methods.  They  are  the  oldest  class  of  solutions,  born  out  of  a  necessity  to  evaluate  early 
astronomical  observations  based  only  on  the  relative  angular  positions  of  wandering  plan¬ 
ets  against  an  unchanging  starfield.  Dr.  Pedro  Ramon  Escobal  provides  an  introduction 
to  such  methods  in  his  1965  text  Methods  of  Orbit  Determination.  He  notes, 

[T]he  angles  only  problem  attracted  the  attention  of  both  Gauss  and  Laplace. 

In  their  day,  this  was  one  of  the  most  pressing  problems  in  mathematical 
astronomy.  Today,  a  century  and  a  half  later,  these  methods  are  widely 
utilized  and,  in  short,  have  stood  the  test  of  time. [Escobal,  1965] 

Should  the  reader  need  a  complete  explanation  of  angles  only  methods  or  common  ce¬ 
lestial  coordinate  systems,  refer  to  Escobal’s  text.  The  rest  of  this  section  explains  only 
the  basics,  because  they  greatly  influence  instrument  design. 

Since  most  orbit  equations  of  motion  use  the  Earth  Centered  Inertial  (ECI)  frame, 
all  measurements  must  be  converted  accordingly.  For  a  fixed  ground  observer,  the  fol¬ 
lowing  parameters  are  required  to  accomplish  this: 

•  Local  Sidereal  Time  (LST):  Since  Babylonian  times,  star’s  longitudes  on  the  celes¬ 
tial  sphere  are  represented  sexigesimally  in  hours,  minutes,  and  seconds  from  an 
arbitrary  point  (the  vernal  equinox).  A  site’s  LST  is  the  celestial  meridian  that 
lies  directly  overhead  at  any  instant  in  time.  To  reinforce  the  fact  this  quantity  is 
an  angle  and  not  time  as  commonly  known,  LST  is  referred  to  as  the  angle  9. 

•  Latitude  and  Longitude:  These  parameters  are  required  in  a  number  of  vector 
transformations,  discussed  in  detail  later. 

•  Altitude:  This  parameter  allows  minor  geometric  corrections  to  the  basic  oblate 
Earth  model  employed  in  transformation  calculations. 

•  Target  Azimuth  and  Elevation:  If  all  previous  parameters  are  available,  a  measure¬ 
ment  in  the  local  horizon,  or  South-East-Zenith  (SEZ)  frame,  may  be  converted  to 
an  inertial  one. 


To  develop  an  initial  orbit  estimate,  position  vectors  are  required.  Gauss  devel¬ 
oped  a  method  of  accomplishing  this  given  only  three  line-of-sight  vectors  and  their 
corresponding  times.2  If  more  than  three  vectors  are  available,  there  are  a  number  of 
algorithms  that  can  improve  solution  accuracy.  Longer  arcs  between  points  are  also 
desirable  [ Escobal ,  1965]. 

In  principle,  angles  only  orbit  determination  is  simple.  Tasks  that  are  simple  in 
theory,  however,  often  become  complex  in  execution.  The  following  section  explores  a 
number  of  complicating  factors. 

2.3  System  Architecture  and  Instrumentation 

Section  2.2  provided  a  short  list  of  data  needs  for  angles  only  orbit  determination. 
This  section  describes  the  system  designed  for  this  project,  assumption  rationale,  and 
other  useful  background  items.  The  architecture  used  in  this  study  is  summarized  in 
Figure  2.3.  Basic  hardware  descriptions  follow,  whereas  specific  calculations  and  software 
components  are  laid  out  in  subsequent  chapters. 
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Figure  2.3:  Angles  Only  Orbit  Determination  Architecture 


2The  “angles”  in  angles  only  refers  to  the  angles  between  line-of-sight  vectors  in  the  orbit  plane.  In 
practice,  the  orbit  is  inclined  from  the  coordinate  frame,  so  an  additional  angle  for  each  measurement 
is  required  in  order  to  determine  the  orbit  plane  orientation.  Therefore,  a  pair  of  angles  (in  this  study 
azimuth  and  elevation)  is  required  to  compute  each  line-of-sight  vector. 
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Before  continuing  to  descriptions  of  specific  hardware  and  software  elements  that 
will  fill  the  architectural  needs  shown  in  Figure  2.3,  an  important  question  must  be 
addressed:  what  precision  and  accuracy  is  required  to  complete  the  task  at  hand?  For 
this  system,  as  with  any  instrument  design,  unbiased  and  random  (i.e.  Gaussian)  errors 
are  desirable.  In  his  book  Modern  Orbit  Determination ,  William  Wiescl  states  that, 

In  practice,  [Gaussian  distribution]  is  achieved  by  finding  and  eliminating 
all  of  the  large  error  sources  in  an  instrument,  until  the  point  of  diminishing 
returns  is  reached.  The  remaining  error  sources  will  be  many  in  number,  and 
small  in  size,  and  the  central  limit  theorem  will  be  obeyed.  [  Wiesel ,  2003] 

Those  familiar  with  experimental  research  know  all  to  well  how  often  equipment  fails  as  a 
white  noise  generator,  however.  To  effect  successful  data  collection,  some  understanding 
of  the  instrument  at  hand  is  required.  This  ensures  suitable  reference  frames  are  selected 
and  appropriate  precautions  taken  throughout  the  design  process. 

Since  this  thesis  examines  commercial  telescopes,  some  aspects  of  the  architecture 
are  predetermined.  For  this  project,  a  homemade  imaging  camera  was  used  in  conjunc¬ 
tion  with  a  popular  commercial  telescope,  the  Meade  LX200GPS  (Figure  2.4). 


10 


Azimuth  Axis 


Handbox 

Direct  Telescope  Control 


Telescope  Computer  and 
Interface  Panel 


Optical  Tube  Assembly  (OTA) 


Coaxial  Wide-Field 
Digital  Camera 


— j—  Elevation  Axis 


Figure  2.4:  Meade  LX200GPS  Telescope 


This  telescope  includes  a  number  of  modern  features,  most  notably: 


•  Global  Positioning  System  (GPS)  Receiver:  Inexpensive  receiver  chips  have  revo¬ 
lutionized  hobby  astronomy,  quickly  calculating  accurate  time,  latitude,  longitude, 
and  altitude  data. 

•  Onboard  Computer:  The  computer  processes  GPS  data,  controls  axis  servomo¬ 
tors,  applies  tip,  tilt,  refraction,  and  other  corrections,  and  has  a  large  internal 
database  of  celestial  objects.  The  telescope’s  computer  communicates  with  a  per¬ 
sonal  computer  via  an  RS232  interface,  using  a  command  set  provided  by  the 
Meade  Corporation  [Meade,  2003]. 

•  Digital  Video  Camera:  A  wide  variety  of  digital  cameras  may  be  used  either  through 
the  telescope’s  main  optics,  or  coaxially  mounted  to  the  telescope’s  Optical  Tube 
Assembly  (OTA).  The  Universal  Serial  Bus  (USB)  architecture  is  commonly  used 
to  both  power  the  sensor  and  transmit  images. 
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In  many  respects,  a  telescope  is  only  as  good  as  its  mount.  This  is  especially 
true  if  an  operator  (or  a  computer)  needs  reliable  attitude  references  in  order  to  find 
targets,  track  them,  or  extract  measurements.  Azimuth-elevation  (or  Az-El)  mounted 
telescopes  like  the  Meade  LX200GPS  are  very  popular  for  all  these  reasons.  Whether 
using  internally-generated  calculations  or  following  ones  from  an  external  source,  the 
telescope’s  onboard  computer  can  independently  command  axes  servomotors  to  recreate 
any  arc  across  the  sky.  For  astronomical  purposes,  this  is  required  to  track  heavenly 
bodies  at  a  sidereal  rate.3  Azimuth-elevation  telescopes  also  suffer  from  a  phenomenon 
known  as  field  rotation.  Objects  appear  to  orbit  around  a  sidereally-tracked  center  point 
in  the  image  plane;  the  effect  gets  worse  further  out.  The  effect  is  noticeable  when 
performing  long-exposure  astrophotography,  but  is  inconsequential  for  exposures  on  the 
order  of  seconds  or  even  minutes. 

Just  as  it  tracks  stars,  an  az-cl  telescope  can  also  follow  satellites  across  the  sky. 
The  Meade  LX200GPS  allows  users  to  upload  satellite  element  sets,  after  which  the 
computer  propagates  their  orbits  and  tracks  them  at  the  appropriate  time.  This  allows 
owners  to  experience  the  “exciting  challenge”  of  satellite  observing  [Meade,  2003].  If  the 
equipment  is  operating  properly,  the  satellite  remains  stationary  in  the  sensor’s  field  of 
view  as  stars  go  whizzing  by  behind  it.  In  theory,  observation  functions  for  this  kind  of 
data  could  produce  accurate  angular  measurements. 

Imagine  the  following,  however:  in  an  American  football  game,  the  ball’s  position 
is  recorded  after  every  play.  Lines  are  painted  on  the  held  to  help  the  referees  estimate 
where  it  stopped,  and  when  necessary,  chains  are  used  to  measure  distances  accurately. 
The  hopes  and  dreams  of  millions  of  fans  sometimes  depend  on  these  measurements. 
Now  imagine  a  game  in  which  the  referees  can  only  watch  the  ball  through  telephoto 
lenses  from  the  top  of  the  stadium,  they  must  keep  the  ball  centered,  and  the  ball  never 
stops  moving.  Placing  the  ball’s  position  this  way  is  hardly  ideal,  just  as  measuring  a 

3 A  sidereal  rate  is  one  revolution  per  sidereal  day,  which  is  «  1/365  shorter  than  a  solar  day.  A 
traditional  equatorial  (or  polar)  mount  aligns  the  telescope  with  constant  lines  of  celestial  latitude,  so  it 
need  rotate  in  only  one  axis  to  keep  a  star  centered  in  the  field  of  view.  Gimballed  mounts  achieve  the 
same  for  satellites,  where  the  telescope  is  first  aligned  with  the  overflight’s  arc,  then  is  panned  either 
left  or  right  at  the  appropriate  angular  rate.  Az-El  telescopes  must  translate  in  both  axes  to  achieve 
the  same  goal,  so  complex  interpolations  are  required. 
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satellite’s  position  against  a  moving  background  is  equally  difficult.  Every  measurement 
requires  a  reference:  for  this  system,  a  relatively  stationary  starfield  will  do  nicely. 

By  astronomical  standards,  this  project  uses  a  relatively  wide-held  digital  camera. 
This  serves  two  purposes:  a)  it  is  more  forgiving  of  targeting  or  timing  errors,  and  b) 
any  given  image  will  be  more  likely  to  have  reference  stars  in  it.  This  is  critical,  because 
mobile  systems  like  this  one  behave  differently  during  every  sortie.  Although  care  is 
taken  to  mitigate  most  major  errors,  it  is  impossible  to  remove  every  misalignment  that 
may  occur.  Since  stars  have  very  well-known  positions,  any  collected  image  may  be 
realigned  using  astrometric  principles.  Details  are  provided  in  Chapter  V. 

The  wide-held  camera  used  in  this  project  is  a  simple  device,  built  from  a  Logitech 
3000  webcam  and  a  vintage  Single  Lens  Rehex  (SLR)  camera  lens.  The  SLR  lens,  a 
35mm  Schneider- Kreuznach  f/2.8,  would  normally  produce  very  wide  (~  70°)  helds  of 
view  when  used  in  conjunction  with  35mm  him.  However,  when  placed  in  front  of  the 
webcam  (with  its  original  lens  removed),  there  is  signiheant  magnification:  the  webcam’s 
chip  is  a  fraction  of  the  size  of  35mm  him.  The  camera  lens  is  positioned  in  front  of  the 
webcam  aperture  so  that  the  system  has  approximate  focus,  and  fine  adjustments  are 
made  using  the  focus  ring  on  the  SLR  lens.  Conveniently,  both  the  webcam  and  lens  ht 
neatly  in  standard  two  inch  PVC  pipe,  which  also  matches  the  guide  scope  mount  on 
the  OTA.  No  great  effort  was  spent  designing  the  sensor,  but  the  following  parameters 
are  empirically  determined  from  using  it  in  practice: 

Table  2.1:  Wide  Field  Camera  Parameters 

Field  of  View  (4:3):  5.7°  /  4.73° 

Image  Size  [pixels]:  640  x  480 
Framerate  [fps] :  5 

Lower  Apparent  Magnitude:  ~  6 

It  is  desirable  to  provide  a  recognizable  starfield  reference,  yet  avoid  complications 
caused  by  a  moving  telescope  such  as  vibrations  or  poor  tracking  performance.  Therefore, 
this  project  adopts  a  slew-and-shoot  method:  a  single  sensor  achieves  hemispherical  sky 
coverage  by  riding  on  a  precisely  aligned  telescope.  The  telescope  leads  the  target, 
pausing  in  anticipation  of  an  intercept.  By  the  time  data  is  collected,  the  scope  is 
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no  longer  moving.  This  method  may  also  be  directly  applied  to  permanently  mounted 
cameras,  either  singly  or  in  a  cluster  configuration.  Regardless  of  the  sensor  setup,  each 
data  collect  (in  this  case  video)  must  have  the  corresponding  metadata: 

•  Time:  Each  video  frame’s  collection  time  enables  precise  dynamic  measurements. 
Since  video  systems  have  relatively  accurate  framerates,  only  the  video’s  start  time 
must  be  logged.  Subsequent  times  are  determined  by  multiplying  the  inverse  of  the 
framerate  by  the  number  of  frames  elapsed  since  video  start.  This  method  assumes 
times  are  logged  in  UTC  time,  then  converted  to  Julian  Dates  using  Equation  3.4. 

•  Site  Location:  For  a  stationary  observer,  site  parameters  must  be  logged  only  once 
per  observing  session.  They  are  used  to  determine  the  site’s  inertial  position  as 
discussed  in  Chapter  III. 

•  Sensor  Altitude  and  Azimuth:  The  telescope’s  reported  attitude  in  terms  of  com¬ 
pass  azimuth  (A)  and  elevation  above  the  horizon  (h)  is  recorded  for  each  video. 

This  document  does  not  include  an  exhaustive  analysis  of  how  observational  errors 
are  dealt  with  in  orbit  solutions  -  refer  to  Wiesel’s  text  for  a  more  complete  analysis. 
Still,  mitigating  such  errors  is  important  and  a  robust  system  must  account  for  them. 
The  following  items  are  are  primary  areas  of  concern,  presented  here  because  they  lurk 
behind  every  calculation  in  this  examination. 

•  Errors  in  Time:  Because  time  is  critical  to  the  transformation  of  observations  from 
local  to  inertial  frames,  access  to  accurate  time  data  is  of  foremost  concern. 

•  Errors  in  Observer  Position:  Without  an  accurate  understanding  of  the  observer’s 
position  on  the  Earth,  local  (SEZ)  to  inertial  (ECI)  transformations  are  not  pos¬ 
sible.  Furthermore,  some  corrections  for  orientation  depend  on  accurate  position 
data. 

•  Errors  in  Telescope  Orientation:  If  an  instrument  is  not  perfectly  aligned  with 
its  assumed  reference  frame,  a  variety  of  errors  may  occur.  For  a  ground-based 
telescope,  there  are  two  basic  but  important  corrections: 
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—  Earth  oblateness:  If  the  Earth  were  a  perfect  sphere,  a  site’s  global  latitude 
would  directly  correspond  to  a  line  of  celestial  latitude.  Due  to  centripetal 
acceleration,  however,  the  Earth  bulges  at  the  equator.  Unless  the  observer  is 
precisely  on  the  equator  or  at  one  of  the  poles,  a  minor  correction  is  necessary. 

—  Tip  and  tilt:  Errors  occur  when  an  a  calculated  local  zenith  vector  does 
not  match  the  true  zenith  vector.  Generally  speaking,  these  occur  when  the 
telescope  mount  is  not  level. 

•  Sensor  Errors:  Whether  looking  through  an  eyepiece  or  using  a  camera,  misalign¬ 
ments  or  rotations  in  optical  systems  induce  additional  errors.  Although  not  an 
error  per  se,  the  effects  of  atmospheric  refraction  must  also  be  accounted  for  when 
interpreting  sensor  data. 

Through  careful  system  design  these  errors  can  be  greatly  diminished,  resulting  in  quality 
observations,  as  shown  in  Chapter  V.  Appendix  A  explains  some  hardware-unique  issues 
encountered  in  this  project. 
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III.  Predicting  Visible  Satellite  Overflights 


Simply  put,  the  sky  is  very  large  and  satellites  are  very  small.  To  make  matters 
worse,  most  telescopes  have  relatively  narrow  fields  of  view.  Colloquially,  this  is 
referred  to  as  looking  through  the  “soda  straw.”  Setting  up  a  telescope  and  waiting  for 
something  to  fly  across  its  field  of  view  would  waste  many  clear  nights.  This  chapter 
describes  how  to  guide  a  sensor  to  appropriate  targets  given  a  catalog  of  satellite  element 
sets  and  brightness  data.  Although  the  ultimate  goal  is  to  generate  orbits  from  measured 
data  without  such  inputs,  solving  this  inverse  problem  first  reveals  many  fundamental 
concepts.  Figure  3.1  depicts  the  elements  that  are  necessary  to  find  and  track  visible 
satellites. 
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Figure  3.1:  Steps  to  Predict  Bright  Satellite  Overflights 


Predicting  when  satellites  will  be  visible  requires  a)  knowing  where  the  satellite  is, 
b)  knowing  where  the  observing  site  is,  and  c)  estimating  how  bright  the  satellite  will 
be.  Sections  3.1  through  3.4  explain  the  methods  and  mathematics  required  to  identify 
bright  satellite  passes.  Section  3.5  showcases  the  integrated  tracking  software  that  was 
developed  during  this  project. 
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3.1  Satellite  Catalogs  and  Orbit  Prediction 

From  this  point  forward,  the  term  satellite  is  used  in  its  purest  sense,  that  is  to  refer 
to  any  body  orbiting  the  Earth.  A  casual  glance  at  the  NORAD  catalog  reveals  most 
tracked  objects  never  were  or  no  longer  are  active  craft.  Even  so,  “what  goes  up  must 
stay  up”  when  dealing  with  orbital  velocities,  so  without  propulsion  or  control  many 
will  remain  in  orbit  for  decades  or  centuries.  In  this  sense,  NORAD  is  the  curator  for 
a  rapidly  spinning  museum,  making  periodic  rounds  to  ensure  each  artifact  is  properly 
labeled.  Generally  speaking,  objects  tracked  during  this  research  project  have  long  since 
passed  from  the  centerfolds  of  trade  magazines.  Discriminating  active  from  inactive 
objects  is  also  well  outside  the  scope  of  this  work.1 

Predicting  overflights  requires  knowledge  of  any  given  satellite’s  inertial  position  at 
any  given  time.  First,  a  reference  estimate  of  the  satellite’s  orbit  at  some  particular  time 
(or  epoch)  is  required.  These  are  called  element  sets,  or  with  NORAD  data  in  particular 
a  Two  Line  Element  Set  (TLE).  Then,  a  means  of  propagating  that  estimate  forward  is 
necessary,  subject  to  appropriate  equations  of  motion  and  perturbations. 

To  use  NORAD  TLE  data,  the  Simplified  General  Perturbations-4  (SGP4)  model 
is  used.  Comprehensive  descriptions  of  the  method  are  available  in  [Vallado  et  al,  2006] 
and  [ Hoots  et  al,  1988],  so  there  is  little  need  to  belabor  the  details  here.  Both  TLE 
data  and  source  code  in  various  formats  is  available  through  Dr.  T.S.  Kelso’s  website 
celestrak.com  [Kelso,  2007]. 

This  project  uses  a  MATLAB®  version  of  the  SGP4  routine  as  adapted  by  Jeff  Beck 
from  code  originally  written  by  David  Vallado.2  It  is  further  modified  to  be  “vectorized,” 
that  is  all  loop  operations  were  replaced  with  matrix  algebra  operations.3  This  improves 
calculation  time  by  about  a  factor  of  16,  allowing  much  better  realtime  overflight  com¬ 
putation. 

1Rapidly  tumbling  objects  exhibit  visible  flash  periods,  so  it  is  possible  to  declare  them  out-of-control 
in  some  cases. 

2MATLAB®  is  a  registered  trademark  of  The  MathWorks,  Inc. 

3The  one  exception  is  the  solution  of  Kepler’s  Equation,  which  is  normally  solved  through  iteration. 
In  this  case,  each  element’s  corrections  are  computed  simultaneously  once  per  loop,  then  individually 
frozen  from  future  updates  as  soon  as  they  reach  convergence.  This  ensures  results  match  the  published 
validation  cases  included  in  [Vallado  et  al. ,  2006]. 
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The  SGP4  algorithm  ultimately  produces  an  inertial  position  and  velocity  vector 
for  any  satellite  in  question.  These  vectors  are  referred  to  as  fsat  and  vsat  in  subsequent 
calculations. 

3.2  Site  Parameters  and  Overflight  Prediction 

Once  a  satellite’s  inertial  position  vectors  are  available,  half  of  the  overflight  prob¬ 
lem  is  solved.  The  second  part  involves  doing  the  same  for  the  ground  site.  Fortunately, 
the  dynamics  are  much  simpler.  This  section  describes  the  required  transformations  that 
conclude  with  the  determination  of  when  and  where  any  given  satellite  will  be  visible  to 
the  observer. 

Of  all  the  required  parameters  in  the  following  transformations,  the  most  important 
is  accurate  time.  This  problem  is  greatly  simplified  by  using  a  GPS  receiver.  Most 
varieties  provide  Universal  Coordinated  Time  (UCT)  or  Zulu  (Z)  time.  Using  any  readily- 
available  formula,  such  as  the  one  included  in  the  Aerospace  Toolbox  for  MATLAB®  , 
UCT  is  converted  to  Julian  date  (JD).  Determining  JD  is  key,  because  first  it  connects 
us  to  satellite  orbit  predictions,  and  secondly  helps  determine  onr  ground  site’s  inertial 
position. 

To  transform  from  the  Earth  Centered  Rotating  (ECR)  frame  to  the  ECI  frame, 
Greenwich  Apparent  Sidereal  Time  (GAST),  expressed  with  the  angle  9g,  is  required. 
This  angle  is  measured  positive  eastward  between  the  vernal  equinox  and  the  Prime 
Meridian,  i.e.  the  location  of  Greenwich,  England.  For  a  given  time,  9g  is  calculated 
using  a  United  States  Naval  Observatory  (USNO)  algorithm.  First,  a  Greenwich  Mean 
Sidereal  Time  (GMST)  is  found  using: 

D  =  ,JD  -  2451545.0 

GMST  =  18.697374558  +  24.06570982441908D  (3.1) 

Then,  a  correction  for  nutation  in  right  ascension  is  applied,  using  a  formula  called  the 
Equation  of  the  Equinoxes  (eqeq).  First,  approximations  for  Mean  Longitude  of  the  Sun 
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(L)  and  Longitude  of  the  Ascending  Node  of  the  Moon  (f^Moon)  are  used  to  determine 
the  Nutation  in  Longitude  (Axj>). 

L  =  280.47  +  0.9856574 

nMoon  =  125.04  -  0.05295474 
A xj)  ~  —  0.000319  sin(nMoon)  —  0.000024  sin(2L) 

After  calculating  obliquity  (e), 

£  =  23.4393  -  0.000000474 

The  correction  term  eqeq  is  found  from 

eqeq  =  Axjj  cos  (e)  (3.2) 

Finally, 

0g  =  GMST  +  eqeq  (3.3) 

For  this  and  any  following  angle  calculations,  the  resulting  answer  is  converted  to 
decimal  degrees  and  wrapped  to  its  appropriate  domain,  in  this  case  [0°,  360°).  According 
to  the  LISNO,  Equation  3.1  loses  one  arcsecond  (~  1.2e~5  degrees)  per  century.  If 
correcting  with  Equation  3.2,  a  maximum  error  of  0.432  arcseconds  and  a  root-mean- 
square  (RMS)  error  of  0.01512  arcseconds  is  expected  [  US  NO,  2007b,  1981].  An  updated 
method  is  available  in  [ Kaplan ,  2005],  but  the  method  presented  here  is  more  than 
sufficient  for  the  selected  application. 

For  a  stationary  observer,  a  site’s  latitude,  longitude,  and  elevation  must  be  de¬ 
termined  only  once  per  observing  session.  Any  GPS  receiver  provides  the  following 
information: 

•  Longitude:  Manufacturer’s  conventions  vary,  so  care  should  be  taken  determin¬ 
ing  which  hemisphere  corresponds  to  a  positive  angle.  For  an  East  longitude  A e 
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(measured  positive  eastward  from  the  Prime  Meridian),  the  site’s  Local  Apparent 
Sidereal  Time  (LAST,  or  6Site)  is  easily  calculated: 

@site  —  9g  +  Ag  (3-4) 

•  Latitude:  Receivers  report  geodetic  (map)  latitude  0,  which  is  used  in  determining 
a  site’s  ECI  position  vector  (Equation  3.6),  the  ECI-to-SEZ  transformation  (Equa¬ 
tion  3.17),  and  the  SEZ-to-A JtKt  transformation  (Equation  4.2)  [ Escobal ,  1965; 
Vallado  and  McClain ,  2007]. 

•  Altitude:  Terrestrial  navigation  references  altitude  above  Mean  Sea  Level  (MSL), 
but  sea  level  and  the  Earth’s  reference  ellipsoid  are  not  coincident.  Most  GPS 
receivers  report  altitude  above  MSL  ( Hmsl )  as  well  as  the  height  of  its  sea  level 
model  above  or  below  the  ellipsoid  (Hgeoid).  The  receiver’s  height  above  the  refer¬ 
ence  ellipsoid,  H ,  is  then: 

H  =  Hmsl  +  Hgeoid  (3.5) 

Once  time  and  navigational  position  are  known,  the  site’s  ECI  position  vector  may 
be  determined  using  Equation  3.6  to  calculate  rsiteA  The  mean  Earth  radius  ae  and 
flattening  parameter  /  are  calculated  from  precise  worldwide  measurements.  GPS  uses 
the  World  Geodetic  System  1984  (WGS-84)  survey  as  a  reference  frame.  Because  /  is 
very  small,  the  WGS-84  standard  provides  its  inverse  [NIMA,  2000].  Two  preliminary 
values  are  calculated:  the  first  is  eccentricity  squared  (e2)  followed  by  the  radius  of 
curvature  in  the  prime  vertical  (A). 


e2  =  2/  -  / 2 


\/l  —  e2sin2((j)) 


4This  is  a  very  common  transformation  from  the  geodetic  to  the  ECR  frame  with  a  positive  counter¬ 
clockwise  rotation  about  the  Earth’s  axis  by  0Site-  The  ECR  transformation  and  other  practical  GPS 
resources  are  referenced  at  [Dana,  2000]. 
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T  site, EC  I 


(  \ 

Y 

yYsite 

Ysite 

>  =  < 

Z site 

(. N  +  H)  cos (9g)  cos (4>)  cos(A e)  —  (N  +  H)  sin(09)  cos (0)  sin(Ae) 
(N  +  H )  sin(#g)  cos(0)  cos(As)  +  (N  +  H)  cos (9g)  cos(0)  sin(Ae) 
(N(  1  —  e2)  +  H)  sin (0) 

(3.6) 


Of  course,  if  a  site’s  parameters  are  well  known,  a  position  vector  may  be  calculated 
at  any  arbitrary  time.  This  is  required  when  generating  future  predictions. 


At  this  point,  both  the  site  and  satellite’s  position  vectors  are  available.  With  this 
data  alone  it  is  possible  to  apply  what  could  be  called  a  binary  brightness  model;  this 
model’s  main  parameter  involves  transmission  losses  due  to  the  amount  of  earth  between 
the  observer  and  satellite.  Plainly  speaking,  simple  checks  are  made  to  determine  if 
the  satellite  is  above  the  local  horizon  or  not.  Let  us  define  a  line-of-sight  vector  Flos 
between  the  site  and  satellite, 

^  LOS  sat  T  site 


Also,  determine  a  unit  vector  that  points  towards  zenith  in  the  ECI  frame  ( z)  as  a 
function  of  (j)  and  0Site: 

' 

COS (4>)  COS (6 site) 

Z  =  {  cos(0)  sin(0site) 
sin(0) 

Finally,  using  the  definition  of  the  dot  product,  the  zenith  angle  z  may  be  calculated. 


z  =  cos 


z-rLOS\ 
\rws\  ) 


(3.7) 


At  this  point,  elevation  from  the  local  horizon  h  is  substituted,  where  h  =  90°  — 
z.  Using  h,  Equation  3.7,  and  either  a  loop  or  vector  of  time,  overflight  occurrences 
are  quickly  calculated.  If  desired,  a  simple  logical  inequality  can  compare  results  to  a 
minimum  threshold,  say  10°  above  the  horizon,  and  disqualify  any  pass  that  fails  to 
break  this  point. 
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This  method  may  also  be  used  to  calculate  sunrise  and  sunset.  The  Jet  Propulsion 
Laboratory  (JPL),  through  their  HORIZONS  online  interface,  provides  information  on 
over  40,000  solar  system  objects.  The  system  can  produce  position  vectors  between  any 
two  bodies  between  selected  times  [JPL,  2008].  Using  Equation  3.7  and  the  Sun’s  ECI 
position  vector  rsun ,  its  zenith  angle  is  determined.  Since  standard  definitions  of  sunrise, 
sunset,  and  twilight  are  only  a  function  of  zenith  angle,  these  times  are  easily  determined 
[USNO,  2007a], 

The  Sun’s  position  vector  serves  a  second  purpose,  as  well.  Although  popular 
media  often  portrays  satellites  as  beeping  behemoths  covered  in  bright  blinking  lights, 
this  is  sadly  not  the  case;  they  emit  no  visible  light  of  their  own.  The  following  method 
determines  whether  or  not  a  satellite  is  directly  illuminated  by  the  sun,  using  a  very 
simple  model  that  treats  the  Earth’s  shadow  as  an  infinitely  long  cylinder  of  Earth’s 
radius.  First,  the  acute  angle  between  the  Sun  and  satellite  vectors,  rjsun,sat  is  found: 


f ]Sun,sat. 


COS 


(  P Sun, ECI  '  P sat, ECI 
\  |  P Sun, EC  1 1 1  Pat, EC  I 


(3.8) 


If  rl Sun, sat  =  0°,  the  satellite  is  directly  between  the  Sun  and  the  Earth.  Conversely, 
rj Sun, Sat  =  180°  indicates  the  Earth  is  directly  between  the  Sun  and  satellite:  the  satellite 
is  in  total  darkness.  The  maximum  angle  at  which  the  satellite  falls  outside  the  Earth’s 
shadow  ( rjmax )  is  a  function  of  the  satellite’s  distance  from  Earth, 


T^-r  +  9°°  (3-9) 

|  T* sat  |  / 

At  any  time  rjsUn,sat  from  Equation  3.8  is  greater  than  the  angle  computed  in 
Equation  3.9,  the  satellite  falls  within  the  Earth’s  shadow. 

By  applying  these  simple,  rapidly-calculated  checks  the  number  of  potential  targets 
is  dramatically  reduced.  Only  satellites  that  have  direct  lines  of  sight  to  both  the  Sun  and 
the  site  in  question  remain,  so  the  basic  overflight  question  is  answered.  These  targets  are 
only  potentially  visible,  however.  The  following  section  describes  how  to  further  improve 
observations  by  generating  an  estimate  of  each  satellites’  brightness. 
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3.3  Satellite  Brightness 

If  a  casual  observer  didn’t  know  artificial  satellites  existed,  spotting  one  in  the 
night  sky  might  prove  difficult  to  explain.  Stars  appear  more  or  less  stationary  over  the 
course  of  minutes,  whereas  meteorites  slash  bright  arcs  in  less  than  a  second.  A  LEO 
satellite,  however,  gradually  appears  out  of  nothingness.  It  may  pulsate,  flash,  or  have 
a  barely  perceptible  tint  of  color.  It  takes  minutes  to  move  methodically  across  the  sky 
in  a  perfect  arc,  then  disappears  as  quietly  as  it  came.  Predictably  determining  when, 
where,  and  with  what  intensity  such  events  occur  is  a  source  of  constant  challenge  for 
satellite  observers.  This  section  describes  how  contemporary  satellite  observers  evaluate 
and  predict  viewing  opportunities. 

Within  contemporary  satellite  observing  circles,  three  individuals  are  widely  asso¬ 
ciated  with  satellite  brightness  predictions:  McCants,  Molczan,  and  Matson.  For  many 
years  now,  their  contributions  continue  to  greatly  aid  the  efforts  of  semi-professional 
satellite  trackers  and  now  their  names  permeate  this  pastime’s  vernacular.  Others  de¬ 
serve  credit  as  well:  there  is  no  intentional  slight  by  failing  to  include  them  here.  The 
composite  method  described  in  this  section  incorporates  multiple  contributions:  it  is 
more  than  sufficient  to  meet  this  project’s  goals.5 

To  predict  satellite  brightness  on  any  given  overflight,  this  project  uses  a  formula 
Robert  Matson  published  (and  often  explained)  online  for  the  benefit  of  the  satellite 
tracking  community  [ Matson ,  2008,  2001].  It  is  presented  in  its  final  form  in  Equation 
3.16  below,  but  some  additional  explanation  is  helpful.  For  now,  consider  this  generic 
logarithmic  equation  which  describes  an  object’s  apparent  magnitude  ( Mapp )  as  a  func¬ 
tion  of  target  distance,  orientation,  and  an  intrinsic  magnitude.  A  satellite’s  intrinsic 
magnitude  is  conceptually  identical  to  a  star’s  absolute  magnitude. 

Mapp  =  Intrinsic  Magnitude  +  Distance  Correction +  Orientation  Correction  (3.10) 

5Higher-fidelity  Iridium  flare  modeling  is  discussed  in  [SeeSat-L  User  Group,  2007].  In  [Henize  et  al, 
1994],  observed  magnitudes  are  compared  to  satellite  Radar  Cross  Section  (RCS)  data.  For  information 
on  the  optical  properties  of  common  spacecraft  materials,  consult  [Culp  and  Gravseth,  1996].  Finally, 
[Kervin  et  al.\  provides  an  overview  of  current  government-funded  research  in  optical  tracking. 
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Generally  speaking,  a  ground  observer  won’t  have  an  estimate  of  satellite  attitude. 
Therefore,  orientation  corrections  can  only  be  a  function  of  Sun-satellite-observer  ge¬ 
ometry.  For  this  reason,  orientation  corrections  are  often  calculated  assuming  satellites 
are  Lambertian-scattering  spheres.  In  this  case,  intensities  of  reflected  sunlight  vary  only 
with  changes  in  phase  angle,  rjphase'- 


hphase 


—  if  f LOS  ‘  Sun  f  sat )  \ 

V I  Flos  I  I  (rsun  ~  rsat )  I  J 


(3.11) 


Phase  angle  varies  from  0°  <  rjphase  <  180°.  When  rjphase  —  180°,  the  object  is  directly 
between  the  observer  and  the  Sun.  Conversely,  when  rjphase  =  0°,  the  object  is  fully  illu¬ 
minated  (assuming  the  Earth  wasn’t  blocking  all  the  sunlight,  of  course).  A  right  phase 
angle  ( rjphase  —  90°)  indicates  a  perpendicular  light  path  between  body  and  observer, 
exactly  the  same  conditions  under  which  a  half-moon  appears. 

Understanding  phase  angle  is  critical,  because  two  major  intrinsic  brightness  cat¬ 
alogs  are  in  common  use;  each  relies  on  a  different  reference  phase  angle.  Accordingly, 
different  formulations  of  orientation  corrections  must  be  applied.  They  are  dubbed  the 
McCants  and  Molczan  Methods,  after  their  creators  and  the  catalogs  that  use  their 
respective  assumptions. 

Michael  McCants’  catalog  of  intrinsic  brightness  is  in  use  since  the  1960’s.  It 
assumes  a  full-phase  reference  angle  ( rjref—  0°)  and  predicts  the  brightest  magnitude 
likely  to  be  observed.  Michael  McCants  explains  he  chose  this  system  “because  I  do 
not  want  to  be  ‘surprised’  that  the  object  is  ‘brighter  than  predicted’”  [ McCants ,  2008b], 
Intrinsic  magnitudes  in  the  McCants  catalog  are  identified  using  Mo°. 

Alternatively,  Ted  Molczan  selected  a  half-phase  ( rjref  =  90°)  definition  for  his 
satellite  catalog.  It  seeks  to  predict  the  likely  average  magnitude  of  an  observed  satellite. 
Many  brightness  values  in  this  catalog  are  based  on  size  and  shape  estimates,  whereas 
others  are  observationally-derived:  the  catalog  annotates  which  method  was  used  for 
each  satellite  [SeeSat-L  User  Group,  Undated].  This  catalog  is  currently  in  the  care  of 
Michael  McCants  [ McCants ,  2008c].  Intrinsic  brightness  values  published  at  Heavens 
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Above  also  use  the  Molczan  Method  [Peat,  2008].  Intrinsic  magnitudes  recorded  using 
the  Molczan  Method  are  identified  using  Mgo°. 

Naturally,  there  is  some  debate  regarding  the  merits  of  each  system.  Robert  Matson 
suggests  the  Molczan  method  is  worthy  for  four  reasons,  paraphrased  here  [ Matson ,  2008]: 

•  Because  T]ref  =  90°  is  in  the  middle  of  the  complete  phase  angle  range,  maximum 
extrapolation  is  only  90°.  Extrapolation  using  the  McCants  method  can  reach  as 
high  as  160°. 

•  For  most  typical  observations,  70°  <  r]p)iase  <  130°,  so  matching  predictions  near 
these  angles  will  produce  better  results. 

•  Should  a  LEO  satellite  actually  be  observed  near  a  full-phase  angle  of  r/phase  —  0°, 
atmospheric  refraction  and  a  lack  of  visible  reference  stars  will  complicate  compar¬ 
isons. 

•  A  true  full-phase  measurement  is  impossible  because  the  satellite  will  be  in  eclipse. 
If  it  could  occur,  there  may  be  significant  boosts  in  brightness  due  to  direct  reflec¬ 
tions  such  as  those  from  solar  arrays. 

To  further  complicate  matters,  there  is  an  approximate  conversion  between  the  two 
systems,  but  Michael  McCants’  own  disclaimers  should  be  consulted  before  proceeding 
[McCants,  2008b]: 

M0 o  =  M90 o  -1.5  (3.12) 

Whether  half-phase  or  full-phase  definitions  are  employed,  a  suitable  orientation  cor¬ 
rection  that  is  only  a  function  of  rjphase  is  still  required.  Robert  Matson  uses  the  term 
phase  factor  to  describe  just  such  a  function  assuming  a  solar-reflecting,  Lambertian 
sphere  [Matson,  2008,  2001].  This  phase  factor  term  is  found  within  square  brackets  in 
Equation  3.13,  where  it  identifies  the  scalar  multiple  increase  or  decrease  in  brightness 
with  respect  to  r]ref  =  90°,  i.e.  the  Molczan  Method.  By  encapsulating  it  in  the  common 
logarithm  conversion  explained  in  Equation  2.1,  it  now  directly  computes  the  expected 
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apparent  magnitude  of  the  satellite.6 


Orientation  Correction^ 0° 


2-5  log10  sin  (rjphase)  +  ( n 


'^'Hpha 

180c 


cos{rjphase)  (3.13) 


Take  note  that  at  rjphase  —  Vref  —  90°,  this  term  equals  zero.  Its  maximum  contribu¬ 
tion  occurs  when  r)p}iase  =  0°,  at  which  an  object’s  apparent  magnitude  would  decrease 
(become  brighter)  by  ~  1.24  due  to  orientation  alone. 

If  the  origin  of  Equation  3.13  is  unclear,  a  similar  result  will  be  re-derived  for  a  full- 
phase  McCants  definition.  It  starts  with  the  equation  for  the  intensity  of  light  scattered 
by  a  Lambertian  sphere  by  a  distant  light  source.7  Irradiance,  reflectance,  and  geometric 
constants  are  lumped,  since  they  will  soon  be  canceled  in  a  ratio  calculation  (noting 
irradiance  is  identical  regardless  of  the  sphere’s  orientation,  provided  it  is  fully  lit).  The 
angle  rj,  at  this  point,  refers  to  any  arbitrary  phase  angle  [Spiro  and  Schlessinger ,  1989]: 


1  = 


constants  sin  {rj)  +  (V 


ny  \ 

180V 


cos  {rj) 


The  ratio  of  these  two  intensities,  one  arbitrary  and  one  with  rjref  =  0°,  becomes: 


/ 


*re/,0° 


constants  sin  {rj)  +  {n  —  cos  {rj)  sin  {rj)  +  {n  —  cosjrj) 

constants  sin  (0°)  +  {n  —  qf^)  cos(0°)  vr 


Then,  after  Equation  3.13,  the  orientation  correction  for  a  McCant  definition  may  be 
expressed. 


Orientation  CorrectioiiQo  =  — 2.51og10 

Now,  this  term  has  zero  contribution  when  r)phase  =  rjref 
magnitude  can  only  diminish  as  expected. 


1  r 


71 


sin  {rjphase)  + 


7T  — 


'ffljpha 


COs{rjphase )J 

(3.14) 


180° 

0°.  As  rjphase  increases, 


6Any  angles  within  trigonometric  functions  are  assumed  to  be  calculated  in  degrees,  to  keep  consis¬ 
tency  with  the  remainder  of  the  paper  (the  original  source  assumes  angles  in  radians).  Scalar  factors 
ranging  from  0  to  tt  are  still  necessary  due  to  the  nature  of  the  phenomena:  7r  has  nothing  to  do  with 
radians  in  this  case. 

7Robert  Matson  adds  that  this  same  equation  holds  true  for  a  Lambertian  cylinder  when  viewed 
broadside.  This  models  rocket  bodies  well,  provided  they  are  not  viewed  endwise. 
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Thankfully,  both  methods  use  the  same  reference  distance  of  1000  km.  Changes 
in  magnitude  due  to  inverse  square  losses  are  expressed  simply,  provided  tlos  is  always 
expressed  in  kilometers  (km): 


Distance  Correctioni000km  =  2.51og10 


\rLos\ 
1000 km 


2 


Slogio  d^LOsl)  -  15  (3.15) 


By  substituting  the  appropriate  corrections  into  Equation  3.10,  a  complete  formula 
for  predicting  satellite  brightness  is  produced.  For  half-phase  Molczan  intrinsic  magni¬ 
tudes,  the  following  equation  as  provided  by  Robert  Matson  is  used  [ Matson ,  2001]: 


Sin  ( Tjphase )  +  (jr  -  ^  )  COS(Vphase) 

(3.16) 

If  using  full-phase  McCants  measurements  the  last  term  in  Equation  3.16  should  be 
substituted  with  the  one  found  in  Equation  3.14,  and  M0°  used  in  place  of  M90 o.  This 
method  appears  valid  because  it  reproduces  results  from  McCants’  Quicksat  program 
almost  identically,  given  the  same  inputs. 

The  McCants  method  was  primarily  employed  throughout  this  project  in  hopes 
of  finding  only  the  brightest  likely  targets.  If  intrinsic  brightnesses  from  the  McCants 
catalog  were  unavailable  for  a  given  object,  values  from  either  the  Molczan  or  Heavens 
Above  catalogs  were  converted  using  Equation  3.12,  then  processed  in  the  exact  same 
manner  as  original  McCants  values.  Although  this  project  did  not  quantitatively  measure 
observed  satellite  brightness,  this  approach  generally  proved  effective. 

3.4  Targeting  in  the  Local  Horizon  Frame 

After  extensive  calculations,  a  short  list  of  potentially  fruitful  opportunities  is  avail¬ 
able.  Chances  are,  the  previous  section’s  calculations  disqualified  a  great  number  of  ob¬ 
jects.  If  our  observing  sensor  (or  vision)  offered  horizon-to-horizon  coverage,  additional 
information  would  be  unnecessary.  Unfortunately  it  does  not,  so  the  following  question 
needs  an  answer:  where  and  when  shall  we  point  our  telescope? 


Mapp  =  M90 o  +  5  log  (\rLOs\)  -  15  -  2.5  log 
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If  previously-mentioned  parameters  are  available,  this  step  is  trivial.  It  requires 
the  satellite’s  position  vector  rsat  and  a  3  x  3  transformation  matrix,  [Csez,eci],  from 
Transformation  14  in  [ Escobal ,  1965].  The  matrix  [C'sez,ect]  is  a  function  of  9site,  and 
<p.  The  SEZ-referenced  vector  is  found  using: 

A 9 at,SEZ  =  [CsEZ,ECl]  ? iat,ECI 


Ssat 

Sx  Sy  Sz 

(  > 

XSat 

Esat 

>  = 

Ex  Ey  Ez 

< 

YSat 

>  (3.17) 

E sat 

EX  Ey  Zz 

Z sat 

\  > 

The  magnitude  of  the  horizon  vector,  \fsatjsEz\,  is  the  slant  range  to  the  satellite. 
The  satellite’s  compass  azimuth  A  and  elevation  above  the  horizon  h  may  be  deter¬ 
mined  from  rsat,SEZ  using  any  quadrant-checking  conversion  from  cartesian  to  spherical 
coordinates.  Venerable  satellite  tracking  programs,  such  as  Quicksat  or  Heavens  Above, 
perform  some  variation  of  the  calculations  presented  thus  far  to  produce  a  tabular  out¬ 
put  of  targets,  times,  azimuths,  and  elevations  of  interest  for  upcoming  sorties  [Me  Cants, 
2008a;  Peat,  2008].  Armed  with  this  data,  the  observer  can  head  out  with  reasonable 
certainty  of  finding  a  satellite. 

3.5  Integrated  Tracking  Software 

Even  if  preparatory  calculations  indicate  dozens  of  targets  will  appear  during  a 
given  evening  or  morning,  any  number  of  obstacles  stand  in  the  way  of  successful  data 
collection.  Foremost  of  these  is  weather.  Even  if  skies  are  partially  clear,  stray  clouds 
can  block  out  portions  of  the  sky.  Other  times,  a  predicted  target  fails  to  appear  at 
the  designated  time  and  place.  Satellites  can  suddenly  wink  out  of  view  if  they  happen 
to  reenter  the  Earth’s  shadow  or  they  fall  into  an  unfavorable  orientation.  For  casual 
observing,  some  of  these  phenomenon  are  quite  entertaining.  When  collecting  data,  they 
are  maddening. 


To  help  overcome  some  of  these  obstacles,  an  interactive  tracking  package  was 
developed  using  MATLAB®  .  At  a  glance,  it  gives  the  observer  a  comprehensive  target  list, 
expected  overflight  paths,  a  realtime  video  feed,  and  the  ability  to  move  the  telescope 
and  log  data  at  will.  If  a  pass  isn’t  working  out,  a  new  target  is  sought  on  the  fly.  Should 
clouds  block  a  portion  of  the  sky,  overflights  in  that  region  are  avoided.  This  human- 
in-the-loop  approach  not  only  prevents  wasted  effort  and  null  collects,  but  is  certainly 
more  entertaining  than  letting  a  computer  pick  every  target. 
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Figure  3.2:  Integrated  Satellite  Tracking  Software 


The  goal  of  this  interface  is  to  allow  the  operator  to  execute  any  part  of  the  tracking 
process  simply  and  clearly.  It  makes  every  calculation  except  brightness  predictions  on 
demand  (i.e.  each  time  a  satellite  is  targeted).  Should  new  data  become  available  during 
a  session,  such  as  an  updated  element  set,  it  could  be  incorporated  immediately.  Other 
features  include: 


•  Target  Selection:  A  strip  chart  in  the  lower  left  features  satellites  of  interest.  It 
regularly  updates,  indicating  what  satellites  are  or  will  soon  be  overhead.  To  aid 
targeting  selections,  graduated  bars  indicate  predicted  brightness.  Once  a  satellite 
is  selected,  a  number  of  target-specific  updates  occur. 

•  Star  Map:  Predicted  or  current  satellite  overflights  are  displayed  in  the  upper  left, 
along  with  the  telescope’s  current  and  targeted  positions.8 

8 Astronomers  may  find  this  map  “backwards,”  because  it  displays  West  on  the  left  rather  than  right 
like  a  typical  star  chart.  It  is  an  indulgence  of  the  terrestrially-minded. 
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•  Video  Display:  The  upper  right  portion  features  a  realtime  video  feed.  A  count¬ 
down  timer  identifies  the  predicted  time  a  satellite  will  cross  the  image  center. 
Clicking  on  the  record  button  logs  the  video  feed  and  writes  metadata  to  a  com¬ 
puter  file. 

•  Telescope  Status:  Just  below  the  video  display,  a  text  readout  shows  the  current 
and  targeted  pointing  angles.  A  simple  algorithm  predicts  the  travel  time  to  move 
between  targets. 

•  Time  and  Position  Data:  Since  the  telescope  is  a  mobile  system,  it  is  helpful  to 
have  a  readout  of  where  the  system  thinks  it  is.  A  GPS  synchronization  feature 
ensures  the  personal  computer  and  telescope  share  navigational  references. 

This  chapter  developed  a  complete  method  for  identifying  bright  satellites  and 
tracking  them.  It  is  wholly  reliant  on  published  orbital  element  sets,  accurate  naviga¬ 
tional  data,  and  either  observational  or  theoretically-determined  brightness  data.  It  is 
also  an  open-loop  process  that  requires  a  human  operator.  In  the  next  chapter,  many  of 
the  transformations  presented  here  will  be  reversed  to  produce  inertial  measurements. 
This  moves  us  one  step  closer  to  develop  our  own  orbit  estimates,  which  in  theory  could 
eliminate  complete  reliance  on  only  the  NORAD  catalog  in  future  tracking  efforts. 
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IV.  Measuring  Satellite  Orbits 


IF  weather,  equipment,  and  a  user’s  basic  competence  combine  to  produce  direct  ob¬ 
servations  of  satellite  orbits,  a  key  challenge  emerges:  how  will  grainy  photographs 
or  videos  of  streaking  satellites  become  accurate  inertial  observations?  Figure  4.1  shows 
the  necessary  elements  required  to  convert  collected  metadata  and  videos  into  the  inertial 
data  necessary  to  perform  initial  orbit  determinations. 


Figure  4.1:  Steps  to  Produce  Inertial  Measurements 

Given  the  tools  in  Chapter  III,  nearly  anyone  could  go  out  at  night,  look  at  a 
predicted  location,  then  find  and  follow  the  tiny  dot  across  the  sky.  If  they  were  asked  to 
describe  when  and  where  they  saw  it,  however,  it  is  unlikely  they  could  produce  useful 
measurements  unless  they  were  specially  prepared  for  the  task  [SeeSat-L  User  Group , 
1998].  The  camera  and  telescope  used  in  this  project  aid  greatly  the  measurement 
process,  because  many  key  parameters  are  automatically  logged.  This  section  describes 
how  measurements  are  made  in  the  camera’s  reference  frame,  corrected  with  reference 
stars  where  available,  and  then  converted  to  inertial  measurements. 
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4-1  Establishing  a  Sensor  Frame  of  Reference 

The  camera  used  in  this  project  is  mounted  co- axially  with  the  telescope,  and  the 
telescope  is  az-el  mounted.  Therefore,  the  Up-Right-Downrange  (URD)  reference  frame 
is  introduced,  which  allows  video  measurements  to  be  converted  back  to  the  SEZ  frame. 
It  is  a  simple  two-axis  rotation  that  requires  knowledge  of  the  sensor’s  azimuth  and 
elevation  in  the  SEZ  frame,  (this  information  is  provided  by  the  telescope’s  computer 
when  logging  data).1 


Figure  4.2:  Up-Right-Downrange  (URD)  Sensor  Reference  Frame 


1  The  World  War  II  anti-aircraft  gun  is  a  good  analogy  here:  the  camera  is  like  the  gunner  always 
peering  through  the  sight.  From  this  point  of  view,  a  target  is  either  up,  down,  left,  or  right  of  the 
barrel. 
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As  evidenced  by  Figure  4.2,  the  sensor  frame’s  horizontal  axis  R  always  remains 
parallel  to  the  horizon  (i.e.  it  is  always  in  the  S-E  plane).  The  U  axis  points  towards 
the  top  of  the  image  plane  (or  “up”  as  viewed  on  a  screen).  The  third  orthogonal  axis 
D  lies  exactly  in  the  image  center,  extending  downrange  of  the  telescope  (colloquially 
known  as  a  boresight). 


With  this  frame  established,  any  unit  vector  in  the  SEZ  frame  is  readily  transformed 
to  the  image  frame: 
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Figure  4.3:  Video  Measurement  Software 


The  SEZ-to-URD  form  is  presented  in  Equation  4.1  because  the  first  use  of  this 
frame  is  to  project  star  locations  onto  collected  videos.  Using  the  methods  described 
in  Chapter  111,  the  angular  positions  of  celestial  objects  such  as  stars  may  be  reduced 
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to  pointing  vectors  in  the  URD  frame  for  any  given  sensor  attitude.2  Then,  a  simple 
Cartesian  to  spherical  conversion  allows  a  determination  of  how  many  degrees  up,  down, 
left,  or  right  of  center  an  object  should  appear  in  any  given  image.3  For  example,  say  the 
frame  transformation  for  a  star’s  inertial  pointing  vector  from  the  ECI  (through  SEZ)  to 
the  URD  frame  indicates  a  star  is  1°  up  and  to  the  right  of  the  reported  image  center. 
This  1°  figure  is  multiplied  by  the  number  of  pixels  per  degree  (a  function  of  camera  held 
of  view  and  resolution)  to  plot  the  star’s  expected  position  directly  on  the  video.  Figure 
4.3  shows  the  measurement  software  developed  to  accomplish  this.  Reference  stars,  with 
their  associated  identifiers,  are  scattered  throughout  the  video  frame.  These  predicted 
star  locations  now  aid  the  satellite  measurement  process,  which  the  next  section  will 
describe  in  greater  detail. 

4-2  Corrections  in  the  Sensor  Frame 

It  is  unlikely  any  collected  image  perfectly  agrees  with  its  associated  metadata,  so 
a  critical  first  step  is  correcting  for  misalignments.  Using  the  software  shown  in  Figure 
4.3,  projected  stars  may  be  lined  up  with  those  collected  on  the  video.  Figure  4.4  shows 
three  common  kinds  of  misalignments:  azimuth  (A^),  elevation  (A/,),  and  camera  twist 

(t)- 

Currently,  the  correction  process  relics  on  a  human  operator  to  identify  necessary 
corrections  in  the  video  measurement  software  in  trial-and-error  fashion.  The  user’s  goal 
is  to  make  all  recorded  stars  appear  within  their  corresponding  theoretically-determined 
projection  (see  Figure  5(b)  for  an  example).  Typically,  the  camera’s  twist  7  is  determined 
first:  comparing  projections  to  any  video  with  more  than  three  stars  usually  produces  a 
reasonable  estimate.4  Then,  the  image  is  translated  left  or  right  (A^),  then  up  or  down 

2This  project  used  formulas  from  [Meeus,  1998]  to  account  for  precession  and  nutation  of  the  inertial 
axis,  as  well  as  proper  motion  and  atmospheric  refraction. 

3Since  pointing  vectors  have  unit  magnitude,  the  vector’s  length  (often  called  p)  is  irrelevant. 

4The  twist  angle  7,  as  determined  from  calibration  videos,  is  used  in  other  observations  where  only 
one  star  is  visible.  Twist  errors  are  static  throughout  a  single  sortie. 
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Figure  4.4:  Misalignments  in  the  Sensor  Plane 

(A h),  until  all  visible  stars  and  projections  match.  Although  it  is  possible  some  alternate 
permutation  of  the  selected  corrections  could  produce  a  match,  this  seems  unlikely.5 

Once  these  steps  are  complete,  a  satellite  streak  is  measured  (see  Figure  5(c))  to 
produce  an  uncorrected  pointing  vector  in  the  uncorrected  URD  frame.  Thanks  to  the 
star  alignment  process,  most  inaccuracies  may  now  be  remedied.  First,  any  camera  twist 
is  canceled  with  a  single-axis  counter-rotation  about  the  image  center  by  the  angle  7. 
Then,  the  azimuth  and  elevation  reported  in  the  metadata  is  modified  by  adding  and 
A h.  This  produces  a  true  apparent  URD  line-of-sight  vector,  which  is  then  transformed 
to  the  SEZ  frame  using  the  inverse  of  Equation  4.1.  Once  in  the  SEZ  frame,  the  effects  of 
generic  atmospheric  refraction  are  subtracted  to  produce  an  airless  azimuth  and  elevation 
measurement.  Now,  the  user  has  the  necessary  information  to  produce  a  single  line-of- 

5Stars’  positions  are  also  updated  frame-by-frame,  allowing  the  user  to  verify  they  are  drifting  at  the 
appropriate  sidereal  rate.  In  this  sense,  this  process  is  three-dimensional  astrometry:  two  dimensions 
are  represented  by  the  apparent  star  position  on  the  image  plane,  and  the  third  dimension  is  time. 
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sight  vector,  pending  one  final  transformation.  The  next  section  describes  this  step  in 
detail. 


4-3  The  Topcentric  Reference  Frame 

Thankfully,  there  is  only  one  more  frame  transformation  standing  between  the 
satellite  tracker  and  an  orbit  determination.  Most  angles  only  methods  make  use  of  the 
Topocentric  Reference  Frame  (ItJtKt),  referred  to  here  as  it  is  in  [Vallado  and  McClain , 
2007].  Simply  explained,  it  is  a  reference  frame  that  is  always  aligned  with  the  ECI 
frame,  but  whose  origin  is  coincident  with  the  observing  site.  Transformation  of  SEZ 
to  It-IfKi  line-of-sight  vectors  is  a  function  of  6site  from  Equation  3.4  and  site  geodetic 
latitude  0. 


cos {6 Site)  cos(90°  -  0)  -  sin(6>Sjte)  cos(0site)sin(9O°  - 
sin(6>Sjte)  cos(90°  -  0)  cos (9site)  sin (9site)  sin(90°  - 

—  sin(90°  —  0)  0  cos(90°  —  0) 


Any  of  the  measurements  produced  using  the  process  described  in  Section  4.2  can 
be  converted  to  topocentric  line-of-sight  vectors  using  Equation  4.2.  Of  these,  the  three 
vectors  with  widest  angular  separation  are  typically  used  to  calculate  orbits:  they  will 
be  referred  to  as  l first,  l mid ,  and  liast,  respectively.  Chapter  V  explains  that  although 
there  are  some  challenges  to  computing  orbits  with  such  measurements,  it  is  possible  to 
get  useful  predictions  given  the  right  conditions. 
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V.  Results  and  Discussion 


THIS  project  began  with  a  simple  goal:  use  a  commercial  telescope  to  track  satellites 
and  determine  their  orbits.  Within  the  span  of  even  a  few  sorties,  the  fundamental 
capabilities  and  limitations  of  a  system  like  this  one  appear: 

•  A  single-site  telescope,  under  the  right  conditions,  can  produce  sensible  initial  orbit 
estimates.  Geometry  dominates  solution  quality,  however.  This  geometry  is  outside 
the  engineer’s  control,  so  ill  effects  must  be  understood  and  tolerated. 

•  This  prototype  system  is  accurate  enough  to  produce  Cartesian  state  vectors  (and 
therefore  classic  orbit  elements)  given  appropriate  observing  geometry. 

•  Investing  in  sensor  precision  does  not  necessarily  yield  matching  dividends  in  a  tele¬ 
scope’s  ability  to  determine  initial  LEO  orbits.  Since  many  LEO  observations  may 
suffer  from  singular  behavior,  using  many  lower-precision  sensors  may  be  better 
than  using  a  few  precise  ones. 

This  chapter  elaborates  on  these  points,  first  by  explaining  observational  limitations 
in  the  angles-only  method.  Section  5.1  explains  this  using  case  studies  from  collected 
data.  Then,  Section  5.2  relates  the  prototype  accuracy  of  this  system  and  demonstrates 
an  initial  orbit  calculation.  Finally,  Section  5.4  shows  how  well  a  case  study  orbit  prop¬ 
agates  to  aid  future  tracking  network  design. 

5.1  Angles- Only  Data  and  the  Great  Circle  Deviation 

It  is  well-established  that  using  only  pointing  vectors  to  determine  satellite  orbits 
is  a  quirky  prospect.  David  Vallado  assesses, 

Gauss’s  method  using  angles-only  data  receives  mixed  reviews  from  the  as- 
trodynamic  community.  The  opinions  range  from  little  concern  because  the 
method  works  best  for  interplanetary  studies,  to  feeling  that  it’s  not  very 
accurate  for  satellite-orbit  determination,  to  reverence  for  the  achievement 
realized  at  a  time  when  data  was  limited.  [  Vallado  and  McClain ,  2007] 

One  key  phenomenon  that  leads  to  enormous  consternation  is  introduced  here,  using  the 
adopted  term  great  circle  deviation.1  Its  effects  certainly  contribute  to  the  ambivalent 

1A  great  circle  is  the  shortest  arc  connecting  two  points  on  a  sphere.  Aircraft  and  ships  use  great 
circle  navigation  to  travel  the  shortest  possible  distance  between  destinations. 
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attitudes  Vallado  notes  above.  It  is  a  counterintuitive  concept  for  most  engineers,  because 
this  phenomenon  dictates  that  regardless  of  instrument  accuracy,  accurate  results  (or 
any  results  at  all)  may  not  be  attainable.  Based  on  data  collected  during  this  project, 
however,  it  is  the  dominant  parameter  affecting  success  or  failure  of  any  attempted  orbit 
calculation  from  a  single  ground  observer.  Therefore,  its  effects  must  be  explored. 

In  Fundamentals  of  Astrodynamics,  the  authors  Bate,  Mueller,  and  White  reference 
work  by  Moulton;  he  found  that  an  angles-only  solution  (in  this  case  the  Laplacian) 
fails  when  “the  three  observations  lie  along  the  arc  of  a  great  circle  as  viewed  from  the 
observation  site  at  [the  middle  time].”  In  an  inertial  sense,  this  occurs  when  the  observing 
site  is  in  the  orbit  plane  during  the  observation  [ Bate  et  al. ,  1971].  Consequentially, 

•  Topocentric  observations  that  form  a  great  circle  will  fail  to  produce  orbit  solutions. 

•  This  condition  is  unstable:  minute  errors  near  this  point  cause  huge  deviations  in 
results.  Solutions  derived  from  near-great-circle  observations  are  suspect. 

•  Deviation  from  the  great  circle  observation  is  proportional  to  the  site  distance  from 
the  orbit  plane. 

Simply  put,  a  telescope  must  have  the  ability  to  look  “down”  on  the  orbit  in  order 
to  observe  its  arc.  Imagine  placing  dozens  of  ants  on  a  table  and  watching  them  wander 
around  on  a  red-and- white  checkered  tablecloth.  It’s  easy  to  trace  the  path  of  any  single 
ant,  as  well  as  determining  if  it’s  on  a  red  or  white  square.  Now,  attempt  the  same  thing 
with  your  eye  at  the  table’s  edge:  it’s  hard  to  tell  the  difference  between  large,  distant 
ants  moving  quickly  or  close,  small  ants  moving  slowly.  Placing  them  on  a  red  or  white 
square  becomes  nearly  impossible. 

This  same  problem  often  occurs  when  observing  LEO  satellites  with  a  telescope. 
Oftentimes,  the  ground  site  is  at  or  near  the  orbit  plane  and  therefore  can’t  see  the 
orbit’s  arc.  As  points  of  light,  the  track  could  be  a  small,  close  ant  or  a  distant  large  one. 
The  orbit  plane  is  clear  (we  know  we’re  looking  at  the  table’s  edge),  but  it’s  impossible 
to  determine  the  satellite’s  range,  or  in  other  words  how  many  red  and  white  squares  are 
between  us  and  the  ant. 
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Keeping  these  phenomena  in  mind,  consider  the  following  case  study:  the  appar¬ 
ent  path  of  two  distinct  satellites  are  carefully  measured  by  a  ground  observer.  An  orbit 
determination  for  each  is  eagerly  computed,  and  the  results  compared  to  those  predicted 
by  the  SGP4  algorithm  derived  from  NORAD  catalog  data.2  For  this  case  study,  re¬ 
sults  from  south-to-north  overflights  of  an  SL-8  rocket  body  (20433)  and  Cosmos  1980 
(19649)  are  presented.  Table  5.1  lists  the  bulk  deviations  from  the  expected  SGP4  pre¬ 
dictions.  The  angles  between  the  observed  and  SGP4-predicted  topocentric  line-of-sight 
measurements  Alfirst,  A lmid,  and  A liast  are  provided  as  an  indication  of  overall  angular 
accuracy  of  the  measurement:  note  the  similar  values.  The  net  error  of  the  calculated 
Gauss/Gibbs  solution  is  computed  by  subtracting  its  computed  position  and  velocity 
vectors  from  their  corresponding  SGP4-predicted  values.  The  magnitude  of  these  errors 
are  listed  as  Er  and  Ev,  respectively. 

Table  5.1:  Observation  and  Orbit  Determination  Errors  for  SL-8  and  Cosmos  1980 

Overflights 

SEZ  Arc  A  l  first  A  Imid  A  liast  hr  (km]  Ev  [m  /  sec] 

SL-8  Rocket  Body  62s  0.127°  0.029°  0.136°  93A  ^892 

Cosmos  1980  67°  0.157°  0.018°  0.161°  0.63  14.82 


Why  is  it  that  one  calculation  produced  useful  results,  but  the  other  was  signifi¬ 
cantly  fouled?  For  a  prototype  system,  initial  guesses  would  involve  timing  or  pointing 
inaccuracies,  miscalculations  of  reference  points,  optical  misalignments,  and  any  other 
number  of  mundane  possibilities.  In  this  case  or  any  other,  attempting  to  correct  time, 
azimuth,  and  elevation  data  or  find  clear  trends  between  them  and  solution  accuracy 
usually  proves  fruitless.  A  clear  answer  emerges  only  after  examining  the  collection 
geometry,  presented  in  Figure  5.1. 


2  Although  the  accuracy  of  NORAD  SGP4  predictions  is  considered  dubious  on  some  scales,  for  this 
example  they  are  sufficient. 
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Observer  Removed  From 
Orbit  Plane  by  0.73°. 


(a)  Overflight  of  SL-8  Rocket  Body  (20433)  on  16  January  2008 


Figure  5.1:  Two  Similar  Satellite  Overflights 


After  inspecting  these  orientations,  it  is  readily  apparent  that  the  SL-8  track  oc¬ 
curred  when  the  angle  between  the  observing  site  and  the  rocket  body’s  orbit  plane  was 
a  scant  0.73°.  This  places  it  very  near  the  table’s  edge.  The  Cosmos  1980  track,  on  the 
other  hand,  enjoyed  over  three  and  a  half  times  more  separation  at  2.73°.  The  calculated 
results  are  no  longer  so  perplexing. 

One  problem  remains,  however.  Identifying  these  very  telling  planar  orientations 
requires  solid  knowledge  of  the  target’s  orbit.  Since  the  goal  is  to  determine  the  orbit 
given  no  other  information,  some  other  method  of  identifying  near-singular  conditions 
is  required.  Moulton’s  original  description  of  the  singular  great  circle  condition  proves 
useful  in  this  endeavor. 

Measurements  are  reduced  to  produce  three  topocentric  pointing  vectors  (l first, 
l mid ,  and  hast)-  Using  only  these  measurements,  a  measure  of  merit  the  author  calls 
great  circle  deviation  (Cgc)  is  proposed;  it  is  computed  in  the  following  manner: 


Find  a  unit  vector  perpendicular  to  the  “great  circle”  plane  that  passes  through 
both  l first  and  hast :3 

hast  X  l first 


lac  — 


hast  X  l first 


Find  the  projection  of  the  middle  topocentric  observation  lmi(i  onto  Ice,  timid'- 


(^mid  (Jmid  '  ^GCj^GC 


•  Find  the  projection  of  lmid  onto  the  great  circle  plane,  bmid. 


1 mid  mid  & mid 


3Formally,  a  great  circle  is  the  circle  formed  by  the  intersection  of  a  sphere  and  a  plane  passing 
directly  through  its  center.  Moulton  referenced  vectors  only  to  the  middle  observation  time,  so  this 

analysis  does  as  well.  The  vectors  lfirst  and  hast >  when  referenced  to  the  site’s  position  at  the  middle 
observation  time,  form  a  great  circle  intersection  with  a  unit  sphere. 
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Compute  the  angle  between  the  middle  observation  and  its  projection  onto  the 
great  circle  plane: 


(cc  =  COS 


-1 


(Lmid  '  braid 


|  O'lnid  | 


Jmid 


(5.1) 


Finally,  it  is  possible  to  determine  why  some  observations  fail  to  produce  useful  ini¬ 
tial  state  vectors  and  others  do,  even  though  the  telescope’s  accuracy  is  constant.  In  the 
course  of  three  sorties,  hundreds  of  measurements  were  made  on  14  distinct  satellites.4 
Figure  5.2  shows  comparative  plots  of  initial  position  and  velocity  estimates  from  these 
observations  ( tiod  and  viod)  versus  ( gc ■  These  observational  results  confirm  theory, 
clearly  identifying  singular  behavior  as  (gc  approaches  zero.  Since  only  a  few  observa¬ 
tions  resulted  in  high  values  of  great  circle  deviation  (they  were  exceptionally  fortunate 
observations),  their  observed  arcs  are  split  up  into  sub-arcs,  and  those  sub-arcs  are  indi¬ 
vidually  analyzed.  This  not  only  fills  out  the  plot,  but  highlights  following:  an  observer 
must  be  outside  the  orbit  plane  to  achieve  high  (gc  values,  but  the  reverse  is  not  true.  A 
highly  separated  observer  catching  only  a  small  arclength  observation  will  produce  mea¬ 
surements  with  low  (gc  values,  so  successful  initial  orbit  determination  is  unlikely.  In 
the  case  of  the  highlighted  Cosmos  1980  and  SL-14  Rocket  Body  (18215)  passes,  results 
would  be  very  poor  if  the  full-length  observation  was  not  collected. 


4See  Appendix  C  for  select  observations  of  these  satellites. 
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0°  0.5°  1°  1.5°  2°  2.5° 


Cgc 

(a)  Absolute  Value  of  Position  Errors  in  True  Range  Direction 


(b)  Absolute  Value  of  Net  Orbital  Velocity  Errors 

Figure  5.2:  Position  and  Velocity  Errors  vs.  Great  Circle  Deviation 
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To  explore  the  impact  of  arclength  on  (qc,  and  add  further  support  to  the  claims 
made  here,  theoretically-determined  azimuth  and  elevation  measurements  computed  as 
in  Chapter  III  are  used  instead  of  collected  ones  in  a  Gauss/Gibbs  angles  only  routine.5 
Again,  using  Cosmos  1980  and  the  SL-14  as  a  case  study,  (gc  is  calculated  as  a  function 
of  the  arclength  of  observation  for  the  same  scenarios  already  discussed.  Figure  5.3  shows 
these  results,  which  confirms  ( gc  is  a  function  of  both  the  arclength  of  observation  and 
separation  from  the  orbit  plane  at  the  middle  observation.  If  it  were  a  function  of 
arclength  only,  the  two  lines  would  be  coincident.  Finally,  it  implies  that  the  closer  to 
the  orbit  plane  an  observer  is,  the  greater  the  observation  arclength  must  be  to  achieve 
the  same  ( gc  value  and  associated  confidence  in  the  initial  orbit  determination. 

100° 

10° 

1° 

o  0.1° 

0.01° 

0.001° 

0.0001° 

0°  20°  40°  60°  80°  100°  120°  140°  160°  180° 

Arclength  of  Observation 

Figure  5.3:  (gc  Versus  Theoretical  Arclength  of  Observation  for  SL-8  Rocket  Body 
(20433)  and  Cosmos  1980  (19649) 

5 When  compared  to  a  line-of-sight  vector  computed  directly  from  fsat  and  rSite,  these  “perfect” 
measurements  produced  a  net  angular  difference  of  0.0005  ±  0.002°.  This  is  probably  due  to  some 
necessary  interpolations.  This  is  about  half  the  angular  accuracy  of  a  Raven  system,  but  sufficiently 
close  to  true  for  this  argument  [ Thrall ,  2005]. 
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The  second  purpose  of  the  case  study  is  to  show  that,  even  with  very  accurate 
measurements,  the  same  basic  phenomena  are  observed.  Figure  5.4  shows  absolute  errors 
in  range  and  velocity  estimates  for  theoretical  data,  just  as  Figure  5.2  showed  them  for 
experimental  data.  The  important  item  to  note  is  the  singular  behavior  near  (cc  —  0°) 
which  confirms  the  experimentally-determined  results.  Future  research  may  explain  why, 
in  these  plots,  accuracy  gets  much  better  shortly  before  it  gets  much  worse.  This  effect 
is  not  observed  in  the  experimental  data  to  date. 

Understanding  the  impact  of  great  circle  deviation  on  a  telescope’s  ability  to  pro¬ 
duce  orbits  is  critical.  It  is  a  sobering  proposition:  both  a  million-dollar  telescope  and 
a  fifty-dollar  camera,  sitting  next  to  each  other  in  the  orbit  plane  of  a  target  satellite, 
are  both  useless  for  determining  state  vectors  by  themselves.  Granted,  a  more  accurate 
instrument  could  produce  better  results  at  lower  values  of  ( qc ,  but  theory  demands  that 
dividends  diminish  exponentially  regardless  of  the  investment. 

It  should  also  be  noted  that  this  effect  is  greatest  for  short  observations,  which 
is  the  norm  for  optical  LEO  tracking.  It  is  exacerbated  by  the  fact  that  these  objects 
are  usually  illuminated  only  when  the  observer  is  close  to  the  orbital  plane  in  the  first 
place.6  When  tracking  higher  altitude  targets,  the  Earth  will  rotate  further  through  the 
target’s  orbital  plane  or  be  further  separated  from  it  in  the  first  place  (as  with  GEO 
satellites,  unless  the  site  is  equatorial).  This  results  in  much  higher  (qc  values  and  fewer 
complications. 


6Extremely  low  elevations  correspond  to  high  orbit  plane  separation.  Except  for  the  highest  latitudes, 
a  low  elevation  observation  in  an  easterly  or  westerly  direction  would  usually  place  the  satellite  behind 
the  Earth’s  shadow  (no  reflected  sunlight)  or  between  the  observer  and  the  Sun  (little  or  no  reflected 
sunlight).  Since  many  LEOs  are  in  highly-inclined  orbits,  the  only  time  these  conditions  are  avoided  is 
when  passing  at  higher  elevations,  which  means  smaller  orbit  plane  separation. 
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(a)  Absolute  Value  of  Position  Errors  in  True  Range  Direction,  Theoretical  Data 


(b)  Absolute  Value  of  Net  Orbital  Velocity  Errors,  Theoretical  Data 


Figure  5.4:  Position  and  Velocity  Errors  vs.  Great  Circle  Deviation,  Theoretical  Data 
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5.2  Prototype  Accuracy  and  Precision 


In  the  previous  section,  the  impact  of  great  circle  deviation  was  emphasized  to 
simplify  the  discussion  of  overall  system  accuracy.  Figure  5.2  shows  a  sharp  cutoff  where 
great  circle  deviation  is  no  longer  the  dominant  source  of  error  («  1°  or  greater  based  on 
experimental  results).  This  cutoff  threshold,  and  the  accuracy  of  any  results  obtained 
beyond  it,  is  now  dominated  by  much  more  linear  angular  and  temporal  parameters. 

As  discussed  in  Chapter  II,  there  are  a  few  major  sources  of  error  in  angles  only 
computations:  time,  site  location,  and  sensor  azimuth  and  elevation.  Accordingly,  a 
telescope  will  have  associated  accuracy  and  precision  for  each.  Accuracy  is  determined 
through  calibration,  whereas  precision  is  dependent  on  measurement  resolution.  This 
section  establishes  these  parameters  for  this  prototype  system.  The  most  important  pa¬ 
rameters  are  azimuth  and  elevation  precision  and  timing  accuracy,  which  are  summarized 
in  Table  5. 2. 7 


Table  5.2:  Angular  and  Timing  Accuracy 

Angular  Accuracy  (3cr):  ±0.05° 

Timing  Accuracy  (3<r):  ±0.5  seconds 


Both  the  timing  and  site  location  accuracy  of  this  system  is  wholly  dependent  on 
the  telescope’s  GPS  receiver.  Individual  chips  used  in  commercial  telescopes  vary,  but 
the  specific  LX200GPS  used  in  this  project  employs  a  Sony  GXB5210.  Due  to  the  fact 
it  reports  time  to  the  nearest  second  and  position  to  the  nearest  second  of  latitude  and 
longitude,  it  is  assumed  to  be  accurate  to  at  least  these  values.8  Once  synchronized, 
timing  precision  is  a  function  of  the  CPU  clock,  software,  and  the  webcam’s  framerate 
output.  Since  negative  effects  of  timing  precision  remain  unobserved,  it  is  considered  to 
have  negligible  impact  on  computed  results. 


7For  reference,  the  original  Moonwatch  teams  first  claimed  to  make  angular  measurements  within 
1°  with  one  second  timing.  With  practice,  some  groups  claimed  six  arcminutes  («  0.004°)  and  0.1 
second  accuracy  [Engle  and  Drummond ,  1965].  Contemporary  telescopes  in  the  Raven  system  achieve 
one  arcsecond  («  1  x  10-5  degrees)  accuracy  [ Thrall ,  2005]. 

8Refer  to  [Sony  Corporation ,  Undated]  for  chip  specifications,  and  Appendix  A  for  more  on  compli¬ 
cations  in  using  the  onboard  receiver. 
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Establishing  angular  measurement  accuracy  is  accomplished  through  astrometric 
correction.  Using  techniques  presented  in  Chapter  III,  inertial  line-of-sight  vectors  for 
any  visible  stars  are  transformed  into  the  URD  plane.  Then,  each  star’s  identifier  and 
location  is  plotted  directly  on  the  video.9  Figure  5.5  shows  an  example  of  a  well-populated 
starheld,  as  well  as  details  of  a  star  and  bright  satellite  as  they  appear  on  a  typical  video. 
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(a)  Well-Populated  Reference  Starheld 
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(b)  Star  With  Correlated  Position 

(c)  Bright  Satellite  Image 

Figure  5.5:  Examples  of  Video  Measurement 


9The  Yale  Bright  Star  catalog  is  used,  which  is  available  at  [ Hoffleit  and  Warren ,  1991]. 
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First  and  foremost,  the  angular  accuracy  of  the  inertial  to  URD  transformations 
and  their  subsequent  display  is  confirmed  by  the  fact  that,  once  minor  corrections  to 
reported  azimuth,  elevation,  and  twist  are  made,  it  is  possible  to  align  all  reference  stars 
with  those  on  the  video.  Additional  confirmation  is  made  by  computing  arclengths  with 
the  software,  then  comparing  results  with  standard  formulas  with  the  same  purpose, 
such  as  those  found  in  [ Meeus ,  1998].  This  accuracy  only  applies  when  reference  stars 
are  visible,  however.10 

In  its  current  configuration,  each  pixel  spans  approximately  0.009°,  so  it  would  be 
tempting  to  use  that  as  a  measure  of  system  precision.  There  is  naturally  some  error 
in  placing  reference  stars  as  well  as  selecting  the  leading  portion  of  a  satellite  streak, 
however.  Since  it  is  likely  that  any  given  user  can  click  inside  or  very  near  the  circle 
as  shown  in  Figure  5(b)  nearly  all  of  the  time,  that  span  should  serve  as  a  suitable  3cr 
precision.  An  angle  of  0.05°  proves  to  be  a  suitable  estimate  based  on  the  span  of  star 
images  and  empirical  analysis. 

5.3  Calculating  an  Initial  Orbit  Determination 

Having  established  both  the  conditions  under  which  initial  orbit  determination  is 
possible  and  the  expected  accuracy  of  the  system,  it  is  possible  to  show  that  this  system 
is  capable  of  producing  a  useful  state  vector.  As  Section  5.1  noted,  one  observation 
produced  a  measurement  with  a  much  higher  £ gc  value  than  any  other  collected  so  far. 
This  observation  corresponds  to  an  SL-14  Rocket  Body  (18215).  If  the  claims  presented 
so  far  are  true,  a  calculated  initial  state  vector  should  be  relatively  accurate.  Figure  5.6 
shows  the  observation  geometry  for  this  overflight. 

Topocentric  measurements  in  the  ItJtKt  frame  are  developed  using  methods  from 
Chapter  IV,  and  a  Gauss/Gibbs  angles-only  routine  from  [Vallado  and  McClain ,  2007] 
is  used  to  determine  the  satellite’s  orbit.  Results  are  computed  in  Monte-Carlo  fashion: 
each  observation’s  nominal  time,  azimuth,  and  elevation  measurements  are  perturbed 
within  the  3cr  values  found  in  Table  5.2.  For  each  perturbation,  new  topocentric  pointing 

10If  no  reference  stars  are  present,  measurements  must  rely  on  mount  accuracy.  Due  to  minor  problems 
as  noted  in  Appendix  A,  only  measurements  with  star  references  are  currently  included  in  this  analysis. 
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Figure  5.6:  Overflight  of  SL-14  Rocket  Body  (18215)  on  3  February  2008 

vectors  l first-,  Imid,  and  liast  are  computed.  For  each  set  of  three  observations,  a  Cartesian 
position  and  velocity  vector  is  generated:  a  mean  solution  with  expected  3<r  deviations 
was  produced  from  these  Cartesian  vectors.  Figure  5.7  shows  a  comparison  between 
the  computed  and  NORAD  orbits.  Table  5.3  compares  state  vector  results  to  those 
extracted  using  the  SGP4  propagator,  followed  by  the  classical  orbital  element  sets  in 
Table  5.3.  When  examining  orbital  elements,  remember  that  large  variations  in  uj  and 
M  are  expected  for  nearly-circular  orbits.  Summing  the  two  produces  the  Argument  of 
Latitude  (/i):  this  parameter  confirms  the  orbits  are  similar  [ Bate  et  al,  1971]. 
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(a)  Angles  Only  Orbit,  Cross-Plane  (b)  SGP4  Orbit,  Cross-Plane 


(c)  Angles  Only  Orbit  at  Q,  Equator  (d)  SGP4  Orbit  at  S±  Equator 
Figure  5.7:  Graphic  Comparison  of  Angles  Only  vs.  SGP4  Orbits 


Table  5.3:  Gauss/Gibbs  Initial  Orbit  Determination  Results  for  SL-14  Rocket  Body 
(18215),  Epoch  3  Feb  08  00:28:32.1Z 


rx 

ry 

rz  [km 

Gauss/Gibbs  Angles  Only 

2784.8  ±  3.66 

4948.4  ±  10.2 

4081.5  ±2.52 

NORAD/SGP4 

2785.6 

4949.9 

4086.6 

vx 

Vy 

vz  [m/sec] 

Gauss/Gibbs  Angles  Only 

1067.0  ±39.8 

4403.5  ±  36.2 

-6047.6  ±  51.4 

NORAD/SGP4 

1076.0 

4384.3 

-6052.0 

Table  5.4:  Orbital  Element  Sets  for  SL-14  Rocket  Body 

Angles  Only  NORAD/SGP4 


Inclination  (i)  82.406°  82.510° 

Right  Ascension  of  the  Ascending  Node  (0)  246.13°  246.06° 

Argument  of  Periapsis  (cu)  104.16°  184.73° 

Eccentricity  (e)  0.00232  0.000955 

Semimajor  Axis  (a)  [km]  7005.41  7002.30 

Mean  Anomaly  at  Epoch  (M)  39.60°  319.26° 

Argument  of  Latitude  (/i)  143.76°  143.99° 
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5.4  Impact  on  Tracking  Network  Design 

Results  shown  here  are  promising.  With  only  a  single  observation,  a  relatively 
accurate  state  vector  is  determined.  This  section  describes  how,  hypothetically,  a  system 
like  this  one  could  contribute  to  a  larger  tracking  network,  or  perhaps  be  duplicated  to 
compose  one  of  its  own.  Although  further  assessments  are  left  to  future  research,  this 
section  examines  how  effectively  the  case  study  orbit  would  have  worked  if  it  put  to 
better  use. 

Any  initial  satellite  orbit  estimate,  even  a  good  one,  will  quickly  lose  utility.  To 
keep  discussions  straightforward,  it  is  common  to  talk  about  an  orbit’s  in-track  and  cross¬ 
track  accuracy.  Of  these,  cross-track  accuracy  is  the  most  important,  since  it  determines 
the  arc  a  sensor  expects  to  intercept  the  satellite  on.  In-track  accuracy  is  generally  less 
important,  because  it  only  affects  how  early  or  late  the  orbit  will  be  along  that  arc.  As 
Chapter  111  described,  most  observers  set  up  an  intercept  position  in  advance  of  the  pass 
and  simply  wait  until  the  target  appears. 

Figure  5.8(a)  shows  the  cross-track  deviation  of  the  computed  orbit  compared  to 
the  reference  solution  as  it  is  propagated  forward  for  nearly  one  day.  Errors  are  greatest 
at  ±90°  from  the  observation  point,  but  pinch  back  together  on  the  opposite  side  of  the 
Earth.  A  more  practical  consideration  is  how  likely  an  optical  sensor  could  reacquire  the 
target,  but  this  is  a  function  of  observer-satellite  geometry  and  the  sensor  itself.  Since 
a  sensor’s  field  of  view  always  plays  a  critical  role,  the  right-hand  axis  in  Figure  5.8(a) 
shows  the  equivalent  angular  separation  of  the  predicted  track  from  the  satellite’s  true 
position,  provided  the  satellite  flew  directly  overhead.  Therefore,  not  knowing  any  other 
parameters,  Figure  5.8(a)  indicates  it  could  take  up  to  a  full  20°  held  of  view  to  guarantee 
catching  the  satellite  within  the  first  day. 

Just  as  cross  track  errors  correspond  in  practice  to  angular  deflections,  in  track 
errors  are  associated  with  the  prediction’s  “lead”  or  “lag”  of  the  actual  satellite  in  time. 
Accordingly,  both  position  errors  and  their  equivalent  timing  errors  are  presented  in 
Figure  5.8(b). 
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These  figures  are  presented  to  aid  additional  analysis  of  telescope  tracking  net¬ 
works.  Although  accuracy  and  precision  of  telescopes  may  vary,  the  fundamental  trends 
presented  here  will  probably  remain  present.  Additional  observations  may  negate  this 
assessment:  in  the  words  of  Norman  Pogson,  “but  this  I  do  not  anticipate.” 
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VI.  Conclusions 


Envelopes  are  important  in  research:  engineers  are  either  expanding  them  or  de¬ 
signing  on  the  back  of  one.  The  two  activities  are  closely  related.  This  paper  is 
much  like  the  proverbial  envelope,  covered  in  descriptions  of  the  basic  pieces  required  to 
assemble  something  much  larger.  It  affirms  many  things  that  are  established  in  theory, 
but  often  unfamiliar  to  most  contemporary  space  professionals.  There  are  other  bene¬ 
fits  to  this  back-to-basics  approach  as  well.  Just  as  great  designs  start  out  as  simple 
scribbles,  elementary  success  in  this  effort  portends  envelope-expanding  capabilities. 

Preliminary  results  answer  the  basic  question  offered  at  the  outset,  namely  “is  it 
possible  to  generate  orbit  predictions  using  commercial  telescopes?”  The  short  answer 
is  yes,  with  the  added  benefit  that  doing  so  requires  no  special  facilities  and  only  a 
fractional  budget  (although  weather  can  be  problematic).  This  is  only  the  first  of  many 
useful  conclusions,  however.  Because  the  equipment  is  simple,  even  nontechnical  indi¬ 
viduals  can  participate  in  research.  Volunteers  with  diverse  backgrounds,  rather  than 
the  author,  usually  operated  the  tracking  software  and  found  it  fun  and  enlightening. 
These  are  exactly  the  emotions  experienced  by  early  Project  Moonwatch  volunteers,  and 
it  is  rewarding  to  see  echoes  of  that  era.  Whether  the  public  ever  contributes  to  satellite 
tracking  with  the  widespread  impact  they  once  did  remains  to  be  seen,  but  it  is  exciting 
to  imagine  the  implications.  Nevertheless,  sincere  hopes  are  offered  that  the  system  de¬ 
veloped  here  will  at  least  help  educate  those  without  space  experience  and  inspire  new 
space-related  research  efforts. 

Such  research,  when  it  occurs,  should  keep  true  to  the  goal  of  seeking  elegant 
solutions  using  simple  tools.  There  is  a  tendency  to  pursue  bigger,  better,  faster,  or 
more  capable  systems  without  pausing  to  consider  why  they  are  necessary.  It  is  also  easy 
to  overemphasize  theory  and  sacrifice  tradecraft  in  the  process.  Working  on  a  system 
like  this  one  provides  a  natural  counterbalance:  if  solutions  become  too  theoretical  in 
nature,  they  are  no  longer  useful.  Likewise,  as  errors  or  data  trends  are  unearthed,  an 
incessant  stream  of  new  theory  is  required  to  explain  them.  If  future  researchers  do  not 
stray  too  far  from  this  locus,  they  can  guarantee  that  any  obtained  results  will  be  of 
practical  use. 
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For  these  reasons,  improving  this  tracking  system  or  expanding  its  capabilities  is 
a  worthy  effort.  Various  advances  are  achievable  within  the  time  and  effort  appropriate 
for  graduate  study.  These  efforts  could  include  the  following: 

•  Orbit /Observer  Analysis:  As  it  stands,  there  is  not  a  complete  analytic  expression 
for  great  circle  deviation  (qc  as  a  function  of  orbit  and  groundsite  parameters.1  Al¬ 
though  the  author  prefers  this  measurement-derived  figure  of  merit,  there  are  other 
methods  that  accomplish  the  same  thing  (refer  to  the  discussion  of  singular  matri¬ 
ces  in  ??).  Regardless  of  the  selected  approach,  it  appears  higher-elevation  passes 
are  less  desirable.  Without  a  more  complete  analysis  of  the  behavior  of  either  it  is 
very  difficult  to  optimize  groundsite  configuration  for  initial  orbit  determination. 

•  Space  Operations:  Inherent  weaknesses  in  single-site  observing  can  be  negated 
with  the  application  of  binocular  tracking,  where  at  least  two  cameras  track  a 
target  to  triangulate  range.  It  is  an  excellent  operations  optimization  problem, 
which  requires  broad  but  basic  modeling  of  tracking  capability  as  a  function  of 
satellite  brightness,  sensor  quality  and  placement,  targeted  satellite  orbits,  weather, 
et  cetera. 

•  Optics  and  Remote  Sensing:  The  camera  used  in  this  project  was  selected  out  of 
convenience.  A  better  camera,  or  more  likely  a  cluster  of  cameras  could  significantly 
improve  LEO  observations.  A  true  optical  “fence”  with  at  least  a  20°  field  of  view 
would  be  a  welcome  development. 

•  Image  Processing:  Automatic  target  identification  techniques  could  render  manual 
postprocessing  obsolete.  This  enables  both  large  collection  volume  and  possibly 
realtime  solutions,  which  has  obvious  benefits  for  SSA  missions. 

•  Systems  Engineering:  Naturally,  the  critical  component  in  any  future  work  is  effec¬ 
tive  systems  engineering.  A  more  detailed  systems-level  analysis  of  the  previously 
mentioned  elements  would  ensure  any  future  systems  meet  user  requirements  and 
are  both  sustainable  and  interoperable. 

^References  to  this  phenomena  are  more  likely  to  be  found  in  archives  that  contemporary  journals. 
Moulton’s  derivation  was  published  in  1914,  and  ground-based  angles  only  techniques  for  LEO  satellites 
are  long  out  of  vogue. 
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It  goes  without  saying  that  any  one  of  these  contributions  provides  an  excellent 
professional  development  opportunity,  especially  for  those  working  in  space-related  fields. 
A  multidisciplinary  team  would  encounter  many  of  the  challenges  faced  by  operators 
today,  albeit  in  a  much  lower  risk  environment  and  with  much  greater  personal  control 
over  system  development.  This  work  proves  exciting  not  because  it  bring  revolutionary 
new  capability.  Instead,  the  revolution  fights  to  reinvigorate  old  concepts  using  effective 
and  inexpensive  new  technology. 
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Appendix  A.  The  Meade  LX200GPS  Telescope 


THis  appendix  presents  a  few  details  unique  to  working  with  the  Meade  LX200GPS 
telescope.  They  are  not  relevant  to  satellite  tracking  in  general,  but  are  presented 
here  to  aid  future  researchers.  Four  main  issues  are  discussed  here: 

•  Complications  using  the  onboard  GPS  receiver:  In  its  factory  condition,  the  tele¬ 
scope  is  not  able  to  provide  entirely  useful  data  to  the  personal  computer. 

•  Discrepancies  in  attitude  reporting:  Although  the  telescope  has  a  reported  whole- 
sky  targeting  accuracy  of  two  arcminutes  (~  0.03°),  there  are  barriers  to  extracting 
attitude  information  with  the  same  quality. 

•  Control  algorithm:  To  accomplish  slcw-and-shoot  tracking,  an  estimate  of  the 
mount’s  speed  is  required. 

•  Imaging  through  the  main  optics:  The  10”  telescope  optics  were  unused  in  this 
project  because,  when  used  with  the  selected  camera,  only  a  1°  held  of  view  could 
be  obtained.  This  was  deemed  too  narrow  for  the  purposes  of  the  project,  but 
some  brief  notes  on  pursuing  this  are  included. 

The  Meade  telescope  used  here  has  a  built-in  GPS  receiver.  Over  the  years,  the 
actual  chip  used  varies:  this  one  uses  a  Sony  GXB5210  in  particular  (discovered  through 
careful  dissection).  The  chip  itself  is  designed  to  output  standard  National  Maritime 
Electronics  Association  (NMEA)  navigation  messages.  However,  the  chip  is  subordinate 
to  the  telescope’s  computer,  so  there  are  a  few  negative  effects: 

•  Message  Filtering:  The  telescope  only  uses  the  GPRMC  message  type  to  syn¬ 
chronize  its  computer.  This  message  effectively  contains  UTC  time,  latitude,  and 
longitude  data  only.  Also,  the  computer  filters  all  other  messages,  so  it  only  pro¬ 
vides  GPRMC  messages  to  the  user  over  the  RS232  port  in  its  factory  condition. 
Richard  Seymour  is  a  prolific  publisher  on  Meade  telescope  firmware,  and  provides 
firmware  patches  that  override  this  limitation  and  fix  many  other  errors.  His  site 
also  contains  more  detailed  explanations  of  GPS  navigation  messages  [ Seymour , 
2008b].  Using  his  4.2g  software  patch,  available  at  [ Seymour ,  2008c],  it  is  possible 
to  gain  access  to  whatever  NMEA  messages  the  receiver  is  willing  to  provide. 
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•  Message  Availability:  Depending  on  the  chip’s  configured  baud  rate,  only  certain 
NMEA  messages  are  broadcast: 

The  GXB5210  can  output  8  different  types  of  sentence:  GPGGA,  GPGLL, 
GPGSA,  GPGSV,  GPRMC,  GPVTG,  GPZDA  and  PSGSA.  If  9600bps 
or  19200bps  or  38400bps  baud  rate  is  set  for  port  setting,  it  outputs 
7  types  of  sentence:  GPGGA,  GPGSA,  GPGSV,  GPRMC,  GPVTG, 
GPZDA,  PSGSA  as  default.  Moreover,  if  4800bps  baud  rate  is  set,  it 
outputs  4  types  of  sentences:  GPGGA,  GPGSA,  GPGSV,  GPRMC  as 
default.  [ Sony  Corporation,  Undated] 

The  telescope  has  the  receiver  chip  hard- wired  to  4800  baud,  so  with  the  4.2g 
patch  the  GPGGA  message  may  be  used  as  well.  Reading  this  message  allows 
altitude  determination.  However,  unless  patch  wires  were  soldered  onto  the  control 
board  and  appropriate  firmware  changes  made,  it  is  not  possible  to  gain  access 
to  the  GPZDA  message,  which  would  provide  UTC  day /month/year  info  as  well 
[Seymour,  2008a].  Therefore,  in  its  current  configuration,  the  user  must  ensure 
the  personal  computer  is  set  to  the  correct  UTC  day,  since  the  telescope  cannot 
provide  this  information. 

•  Lag  and  Precision  Issues:  Probably  due  to  a  lag  in  the  telescope’s  computer,  GPS 
messages  are  not  streamed  sufficiently  fast  to  be  considered  realtime.  Comparisons 
of  loghles’  mean  determinations  with  their  deviations  show  errors  as  high  as  1.5  ± 
0.009  seconds.  A  working  correction  of  1.5  seconds  is  now  added  to  synchronization 
routines,  but  if  greater  timing  precision  is  required  it  is  recommended  to  abandon 
use  of  the  onboard  receiver. 

In  order  to  target  the  camera  and  record  critical  metadata  for  the  measurement 
process,  the  telescope’s  computer  must  be  able  to  slew  to  and  accurately  record  any 
position  from  horizon  to  horizon.  The  computer  is  optimized  to  handle  only  the  former. 
It  takes  a  desired  target,  usually  a  Right  Ascension  (RA,  a)  and  Declination  (Dec,  5) 
value  from  a  celestial  object  catalog,  and  corrects  for  measured  tip,  tilt,  precession, 
nutation,  and  refraction  effects. 

In  early  attempts  at  pointing  the  telescope,  GoTo  Az  and  GoTo  El  targets  were 
commanded  (:  Sz#  and  :  Sa #  in  the  Meade  Command  Protocol).  The  inverse  measure- 
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merits  (:  GZ #  and  :  SA #)  were  recorded  with  videos.  This  corresponds  to  the  terrestrial 
tracking  mode,  where  the  computer  does  not  automatically  engage  sidereal  tracking  af¬ 
ter  each  slew  is  complete.  Apparently,  it  does  not  account  for  tip  or  tilt,  either.  To 
counter  this,  an  attempt  was  made  at  back-calculating  azimuth  and  elevation  from  the 
telescope’s  reported  RA,  Dec,  computed  LAST,  and  navigational  data  using  transfor¬ 
mations  from  [ Escobal ,  1965].  When  working  in  RA  and  Dec  coordinates,  the  telescope 
assumes  celestial  tracking  mode,  and  automatically  applies  a  number  of  corrections,  as 
well  as  automatically  engaging  sidereal  tracking  after  every  slew. 

In  the  following  experiment  comparing  the  az-el  versus  R A/Dec  methods,  inter¬ 
esting  trends  are  noted.  The  telescope  was  carefully  leveled  (a  carpenter’s  level  was 
placed  on  the  OTA  through  full  travel)  and  slewed  in  a  complete  arc  at  0°,  45°,  and  85°, 
respectively.  Figure  A.l  shows  the  results. 
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Figure  A.l:  Back-Calculated  and  Directly  Reported  Azimuth  and  Elevation 
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These  trends  imply  a  few  things.  First,  the  scope  may  be  subtracting  refraction 
effects  in  its  output  position  using  the  RA/Dec  method.1  The  sine  trend  could  be  due 
to  multiple  things: 

•  LAST  Miscalculation:  Many  observers  note  their  telescope’s  computation  of  LAST 
differs  from  what  it  should  be  (Richard  Seymour  confirms  there  are  both  past  and 
current  examples  of  this).  This  would  result  in  a  net  “twist”  in  results. 

•  Invisible  Precession/Nutation/Tip-Tilt  Parameters:  Richard  Seymour  also  con¬ 
firms  that  the  computer  uses  its  own  internal  coordinate  reference  frame,  which  is 
generally  unavailable  to  the  user.  It  is  possible  this  twist  is  due  to  rotations  caused 
by  “invisible”  parameters  such  as  precession,  tip,  and  tilt.2 

These  issues  may  not  need  resolution,  however,  because  the  RA/Dec  method  turns  out 
to  be  much  slower  than  the  az-cl  one.  Furthermore,  it  is  harder  to  predict  how  long 
a  slew  will  take,  because  the  telescope  overshoots  the  target,  then  approaches  it  in  the 
sidereal  tracking  direction  to  avoid  gear  backlash  issues.  This  is  excellent  for  observing 
stars,  but  detrimental  to  satellite  tracking,  since  it  results  in  lost  observing  time. 

Speaking  of  observing  time,  the  slew-and-shoot  technique  requires  an  estimate  of 
how  long  it  takes  the  telescope  to  move  from  any  arbitrary  point  to  any  other.  Lab 
tests  revealed  that,  when  powered  by  an  external  supply,  a  commanded  slew  moves  at  a 
fairly  linear  4.7—  in  both  axes.  This  is  preceded  by  a  two  second  command  delay,  and 
roughly  a  six  second  settling  time.  These  estimates  are  working  well  in  practice,  since 
the  telescope  will  reach  its  intended  target  (usually  a  five  second  lead  in  the  satellite 
track)  and  settles  out  consistently.  The  one  issue  not  yet  addressed  is  the  occasional 
“unwinding”  that  occurs  when  the  telescope  computer  decides  to  move  its  azimuth  axis 
closer  to  home  position  (this  axis  has  ±360°  net  travel).  A  running  tally  of  slews  could 
be  kept  in  the  software  to  predict  when  such  a  move  is  likely. 

To  date,  the  original  az-el  method  is  being  used.  By  leveling  the  telescope  a  little 
more  carefully,  and  using  only  measurements  with  reference  stars  in  them,  reliable  re- 

1At  the  horizon,  exo-atmospheric  objects  appear  «  0.5°  higher  than  they  actually  are. 

2It  is  possible  failure  to  correctly  apply  the  equation  of  the  equinoxes  is  at  fault,  a  fact  discovered 
shortly  before  publication. 
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suits  may  be  obtained.  Subsequent  researchers  will  undoubtedly  resolve  any  outstanding 


issues. 

Finally,  an  all-too-common  question  is:  “why  do  yon  have  that  big  telescope  if 
yon  don’t  use  it?”  The  answer  is  simple:  a)  buying  the  mount  and  telescope  separately 
isn’t  cost  effective,  b)  the  main  optics  still  enable  precise  alignment,  and  c)  it  may  be 
used  in  later  projects.  The  current  camera  was  tested  with  a  Meade  f/3.3  focal  reducer 
to  determine  what  kind  of  images  could  be  expected.3  Image  quality  was  similar  to 
those  produced  with  the  SLR  lens,  but  the  field  of  view  was  on  the  order  of  1°.  It  was 
deemed  unsuitable  for  this  first  research  attempt,  but  is  certainly  worth  employing  for 
appropriate  projects,  such  as  satellite  tracking  in  higher  orbits.  Better  cameras  should 
definitely  be  employed  whenever  possible,  as  well. 


3A  focal  reducer  is  a  specially-crafted  lens  that  effectively  shortens  the  focal  length  of  the  telescope, 
which  gives  it  a  wider  field  of  view  (to  offset  magnification  caused  by  small  CCDs)  and  focuses  more 
light  on  each  pixel,  resulting  in  faster  exposures. 
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Appendix  B.  MATLA 5®  Function  Descriptions 


NUMEROUS  MATLAB®  scripts  and  functions  are  used  in  the  course  of  this  research 
project.  This  appendix  is  divided  into  four  parts:  one  each  for  the  target  identifi¬ 
cation  script,  integrated  tracking  software,  measurement  software,  and  angles  only  orbit 
determination  tools.  Since  there  are  numerous  inputs  and  outputs,  most  of  which  are 
graphics  handles,  these  are  omitted  for  clarity. 

B.l  Target  Identification  Script 

This  script  performs  all  calculations  presented  in  Chapter  III.  A  user  must  provide 
the  following  input  hies: 

•  NORAD  Three-Line  Element  Set:  The  three-line  version  of  a  NORAD  ele¬ 
ment  set  includes  the  satellite’s  common  name  on  the  line  preceding  the  two-line 
data  set.  Operators  always  want  to  know  what  they’re  tracking.  These  are  avail¬ 
able  from  celestrak.com,  with  filename  format  catalog  SI  JYYYY  _M  M  _D  D  Mm.txt 
or  catalog  Ail  APYYY  _MM  _DD_pm.txt,  depending  if  it  is  the  morning  or  evening 
(Colorado  time)  catalog  release. 

•  qs.mag:  This  is  Mike  McCants’  intrinsic  satellite  brightness  catalog,  posted  at 
[. McCants ,  2008c]. 

•  mcnames:  The  original  Ted  Molczan  catalog,  available  at  the  same  location  just 
referenced. 

•  ha.txt:  An  intrinsic  brightness  hie  generated  from  the  satellite  database  at  Heavens 
Above  [Peat,  2008]. 

•  sun.txt:  Solar  ECI  vectors  from  JPL’s  HORIZONS  website  [JPL,  2008].  Data 
must  be  extracted  with  the  following  options: 

—  Ephemeris  Type:  VECTORS 

—  Target  Body:  Sun  [Sol]  [10]  (note  planets  could  be  added  to  the  current  star 
chart  substituting  planetary  data  here) 

-  Coordinate  Origin:  GEOCENTRIC  (500)  (i.e.  ECI) 
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—  Time  Span:  Choose  any  span  (current  file  valid  through  2008),  but  TIME 
STEP  =  Id  (starting  at  UTC  00:00:00). 

—  Table  Settings:  Output  Units  =  KMSzKM/d 

—  Quantities  Code=  2  (State  Vector:  x,y,z,vx,vy,vz) 

—  Reference  Units:  Earth  mean  equator  and  equinox  of  reference  epoch. 

-  Reference  System:  ICRF/J2000.0 

—  Correction:  None 

—  Labels:  Disabled 

-  CSV  Format:  Disabled 

—  Display/Output:  Download/Save  (this  generates  an  output  .txt  hie). 

Using  this  or  equivalent  source  data,  the  following  output  hies  are  generated: 

•  precalc_results.mat:  This  contains  results  of  brightness  calculations.  Its  most 
notable  variable  is  obs  keep ,  which  lists  satellite  brightness  as  a  function  of  JD. 
Only  satellites  that  beat  a  minimum  threshold  at  one  point  during  the  selected 
timespan  (usually  sunset  to  sunrise).  Dim  passes  are  hltered  by  the  integrated 
tracking  software.  Any  time  the  satellite  is  below  the  horizon,  its  brightness  is 
listed  as  “99.” 

•  precalc_tle.txt:  This  hie  is  unused  by  later  scripts,  but  is  a  trimmed  three-line 
catalog  containing  only  a)  satellites  with  periods  less  than  some  desired  number 
(225  minutes  is  the  default  for  LEO)  and  b)  only  satellites  with  brightness  data  in 
one  of  the  three  catalogs  listed  above. 

•  QUICKSAT.DAT:  This  is  a  Quicksat-friendly  element  set  hie  including  only 
bright  satellite  ephcmerides.  Quicksat  uses  non-Windows  carriage  returns  and  has 
a  maximum  catalog  input  of  3000  satellites  [ McCants ,  2008a]. 
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There  are  many  scripts  used  to  process  these  inputs  and  save  outputs,  most  of  which 
are  inherent  to  core  MATLAB®  code  or  its  Aerospace  or  Mapping  Toolboxes.  Notable  ones 
include: 

•  Vectorized  SGP4  Scripts:  The  following  are  vectorized  versions  of  Jeff  Beck’s 
MATLAB®  adaptation  of  David  Vallado’s  SGP4  code  written  in  C.  The  complete 
package  may  be  found  at  [Kelso,  2007].  Note  only  the  following  hies  are  modified, 
but  they  may  still  depend  on  subfunctions  included  in  the  referenced  collection. 
These  scripts  in  concert  produce  fsat. 

—  dpper .vectorized. m:  Modified  for  vector  output. 

—  dspace.vectorized.m:  Modified  for  vector  output. 

—  sgp4_vectorized.m:  Core  SGP4  code. 

—  sgp4_init_vectorized.m:  Minor  change,  now  points  to  .vectorized  versions 
of  code  so  originals  may  be  deleted. 

—  twoline2rv_simple.m:  Corrects  a  number  of  errors  and  produces  initial  or¬ 
bital  elements  and  epoch  times  for  a  single  satellite  at  a  time,  given  two-line 
input  strings  from  the  catalog. 

•  getcatalogsats.m:  Using  input  data,  produces  precalc_tle.txt. 

•  getdarkness.m:  Computes  sunrise  and  sunset  times. 

•  getstars.m:  Imports  data  from  the  Yale  Bright  Star  catalog  and  produces  a  table 
of  star  identifier,  right  ascension,  declination,  and  apparent  magnitude. 

•  getsite.m:  Calculates  the  site’s  EC1  position  vector,  rs*te- 

•  getsun.m:  Extracts  ECI  sun  vectors  from  sun.txt,  then  interpolates  them  as  nec¬ 
essary  to  get  fSun- 

•  getzenith.m:  Calculates  the  site’s  zenith  vector  in  ECI  coordinates. 

B.2  Integrated  Tracking  Software 

The  integrated  tracking  software  is  shown  in  Figure  3.2.  By  default,  it  refreshes  all 
data  at  a  set  rate  (currently  15  seconds)  unless  otherwise  noted.  For  certain  operations, 
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it  automatically  places  important  information  in  a  file  named  logfile.txt.  This  Graphic 
User  Interface  (GUI)  has  the  following  features: 

•  Comprehensive  Starmap:  The  upper-left  corner  features  a  starmap  that  shows  not 
only  the  current  position  of  numerous  bright  stars,  but  also  target  satellite  paths, 
the  telescope’s  current  position,  and  the  telescope’s  intercept  position. 

•  Brightness  Predictions  and  Satellite  Targeting:  Results  stored  in  precalc -results .mat 
may  be  viewed  graphically  in  the  screen’s  lower  left  corner.  For  any  given  minute 
(the  default  precalculation  step  size),  a  satellite’s  expected  brightness  is  indicated 
with  one  of  six  bar  intensities,  one  each  for  apparent  magnitudes  above  six.  Cur¬ 
rently,  satellites  are  sorted  so  those  with  the  highest  mean  brightness  in  the  next 
15  minutes  are  at  the  top.  By  clicking  on  any  satellite,  it  is  automatically  targeted. 
When  targeted,  a)  the  satellite  path  is  shown  on  the  starmap  in  the  upper  left 
corner,  b)  a  telescope  slew  is  generated  along  with  an  estimated  time  of  travel,  and 
c)  a  countdown  to  the  time  the  satellite  will  cross  the  image  center  begins.  Slots 
are  reserved  for  any  satellite  becoming  visible  in  the  next  25  minutes,  but  they  are 
only  shown  if  rising  in  the  next  15  minutes.1  If  a  satellite  is  not  yet  risen,  an  arrow 
indicates  the  direction  of  its  pass.  Once  a  satellite  is  overhead,  its  future  path  is 
shown  in  red,  its  current  position  is  marked  with  a  square,  and  its  location  within 
the  last  five  minutes  is  shown  in  gray. 

•  Guide  Star  Targeting  and  Identification:  The  “Stars”  button  brings  up  a  menu  of 
currently  visible  guide  stars  (stars  that  the  Meade  telescope  uses  in  its  alignment 
process).  Selecting  one  initiates  the  same  process  as  with  satellites,  except  the 
star  changes  color  on  the  starmap  instead  of  producing  an  overflight  path.  This 
feature  helps  non-astronomers  find  requested  stars  when  asked  for  them  by  the 
Meade  computer.  It  is  also  useful  for  checking  telescope  alignment  and  recording 
calibration  videos. 

•  Video  Preview  and  Recording:  The  upper  right-hand  corner  shows  a  preview  of  the 
camera  video,  with  a  central  crosshair  superimposed  over  it.  Immediately  below 

^Reserving  slots  prevents  the  chart  from  becoming  to  “jumpy”  as  new  sats  appear.  Otherwise,  a 
satellite’s  ranking  may  change  rapidly. 
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it,  a  box  displays  the  UTC  time  a  satellite  is  expected  to  cross  the  image  center 
(assuming  the  telescope  reaches  intercept  by  its  predicted  time),  a  countdown  to 
the  same  time,  and  a  “Record”  button  that  logs  both  streaming  video  and  relevant 
metadata  to  the  loghle.2 

•  Telescope  Slew  Controls  and  Data  Readout:  A  text  output  of  selected  target,  esti¬ 
mated  time  to  intercept,  current  telescope  position,  and  intercept  telescope  position 
is  shown  below  the  video  recording  box.  The  Slew  command  button  initiates  an 
immediate  recalculation  of  all  orbit  and  intercept  geometries,  commands  the  tele¬ 
scope  to  slew,  and  automatically  deselects  the  button  when  a  slew  is  complete 
(once  the  telescope  reports  the  same  position  for  three  consecutive  seconds).  The 
Track  and  Sync  buttons  are  mutually  exclusive. 

•  Location  Information:  As  a  check,  location  information  is  provided  in  the  lower- 
right  hand  corner.  The  Sync  button  initiates  a  synchronization  routine  with  the 
telescope’s  onboard  computer.  It  also  logs  raw  GPS  message  data  to  the  loghle. 
The  Sync  and  Track  buttons  are  mutually  exclusive. 

•  Time  Selection  Box:  Normally,  the  software  is  used  in  realtime,  which  requires 
selecting  the  Current  Time  radio  button.  To  view  any  other  time,  the  operator 
may  select  User  Selected  Time,  which  allows  examination  of  scenarios  in  the  past 
or  future.  Typically,  this  feature  is  used  to  look  ahead  during  a  sortie  to  find  out 
if  suitable  targets  still  remain. 

As  before,  numerous  functions  are  required  to  generate  the  GUI  and  process  user 
inputs.  Significant  scripts  include: 

•  trackgui.m:  Master  script,  generates  the  majority  of  interface  objects  and  execu¬ 
tion  timers. 

•  refreshgui.m:  Triggers  GUI  update  on  a  set  refresh  cycle. 

•  createschedule.m:  Generates  the  brightness  data  strip  chart  and  triggers  target- 
related  updates  when  a  satellite  is  selected  from  the  plot. 

2The  “Auto”  label  refers  to  an  unimplemented  feature  that  automatically  begins  logging  video  when 
the  countdown  timer  reaches  a  set  threshold. 


68 


•  createstarmap.m:  Computes  expected  positions  of  visible  stars  and  plots  them 
on  the  starmap. 

•  buildsattrack.m: Recalculates  the  expected  time,  azimuth,  and  elevation  of  a 
satellite  overflight  and  plots  its  path  on  the  starmap. 

•  findtelescope.m:  Queries  available  ports  for  an  active  Meade  telescope,  then 
automatically  connects  if  one  is  found. 

•  getscopeposition.m:  Queries  the  telescope  for  its  current  position.  If  errors  are 
encountered,  it  generates  a  “No  Scope”  error  in  the  tracking  software. 

•  gpssync.m:  Imports  the  telescope’s  GPS  data  stream  and  logs  relevant  naviga¬ 
tional  data.  It  runs  a  comparison  of  the  reported  stream  time  with  the  computer’s 
own  clock,  then  notes  the  offset.  If  the  standard  deviation  is  sufficiently  low,  the 
tracking  software  applies  this  offset  to  all  subsequent  time  calculations  and  dis¬ 
plays.  If  not,  a  warning  occurs.  If  other  errors  are  encountered,  it  generates  a  “No 
Scope”  error  in  the  tracking  software. 

•  slew.mcancelslew.m:  Initiates  or  cancels  telescope  slewing  while  employing  a 
number  of  crosschecks.  If  errors  are  encountered,  each  generates  a  “No  Scope” 
error  in  the  tracking  software. 

•  slewbuilder.m:  Given  a  current  telescope  location  and  a  targeted  satellite,  this 
function  estimates  the  approximate  travel  time  required.  It  then  returns  a  target 
intercept  position  and  estimated  time  the  satellite  will  cross  image  center. 

B.3  Video  Measurement  Software 

Video  measurement  software,  shown  in  Figure  4.3,  allows  a  user  to  import  collected 
videos,  correct  for  misalignments,  and  measure  satellite  streaks.  It  automatically  reads 
loghle.txt,  generated  with  the  integrated  tracking  software  during  observations.  Once  it 
imports  the  metadata  associated  with  each  video  (and  if  every  piece  of  critical  metadata 
is  present),  the  software  calculates  the  expected  position  of  any  visible  stars  and  displays 
them  directly  on  the  video.  A  Show  Stars  checkbox  removes  computer-predicted  guide 
star  locations.  Using  controls  on  the  right,  corrections  for  azimuth,  elevation,  and  camera 
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twist  are  interactively  applied.  Once  this  step  is  complete,  a  user  need  only  click  on 
the  satellite  streak  to  generate  a  JD  time,  azimuth,  and  elevation  record  under  the 
“measurement”  variable. 

Major  scripts  and  functions  include: 

•  measurevideo.m:  Main  script  that  generates  the  GUI  and  defines  all  button 
callback  functions. 

•  getstars.m:  As  previously  described.  Assumed  epoch  is  the  computer  date/time, 
stars  are  trimmed  to  apparent  magnitude  six  and  above  by  default.  This  function 
runs  only  once  at  script  startup. 

•  loadvideo.m:  Imports  video  data,  erases  bad  pixels  by  matching  their  intensities 
to  the  frame  mean,  and  identifies  corresponding  metadata  from  the  loghle. 

•  getstarazel.m:  Converts  stars’  ECI  pointing  vectors  from  getstars.m  to  the  SEZ 
frame  using  each  video’s  metadata,  then  trims  all  results  to  only  those  stars  above 
the  horizon  after  accounting  for  generic  atmospheric  refraction.  Repeats  on  each 
frame  refresh,  allowing  the  stars  to  drift  at  their  appropriate  rates  (assuming  a 
stationary  camera).  Currently,  every  star  is  transformed  in  this  step,  because  no 
serious  delays  are  observed. 

•  buildfovstars.m:  Further  trims  stars  to  include  only  those  that  will  appear  in 
the  URD  frame  (i.e.  the  video),  determines  their  corresponding  locations  in  pixels, 
then  plots  and  labels  them  on  the  video. 

•  measureimage.m:  Whenever  the  video  is  clicked,  this  function  identifies  the  pixel 
clicked,  converts  it  to  a  corresponding  URD  pointing  vector,  corrects  it  using  the 
azimuth,  elevation,  and  twist  errors  identified  earlier,  transforms  it  to  the  SEZ 
frame,  then  unrefracts  it  to  an  airless  measurement. 

B.4  Angles  Only  Orbit  Determination  Tools 

Once  airless  JD  time,  azimuth,  and  elevation  measurements  are  available  from  the 
measurement  software,  they  may  be  used  in  conjunction  with  site  parameters  to  produce 
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topocentric  pointing  vectors  and  subsequent  orbit  solutions.  A  few  variations  of  scripts 
are  used  to  not  only  determine  initial  orbits,  but  also  analyze  and  interpret  results. 
This  section  describes  the  baseline  deterministic  configuration:  measurements  are  used 
directly,  so  no  uncertainty  is  generated.  The  entire  method  follows  from  [  Vallado  mid  Mc¬ 
Clain ,  2007],  which  should  be  consulted  for  details.  This  script  measureconverter.m 
requires  the  following  inputs: 

•  measurement:  This  variable  is  produced  using  the  video  measurement  software. 
It  is  typically  imported  from  NNNNN.mat,  a  workspace  file  containing  the  mea¬ 
surement  variable  for  a  given  satellite  ID. 

•  azelrange:  This  variable  is  the  result  of  using  the  satpath.m  script,  which  is  a 
version  of  precalcs. m  that  computes  observation  geometry  for  only  one  satellite.  It 
is  a  table  of  JD  times  with  a  corresponding  azimuth  and  elevation  for  each.  This 
variable  allows  error  calculations  and  a  display  of  “true”  results  for  comparisons. 
Otherwise,  it  is  unnecessary. 

Given  these  inputs,  the  script  automatically  selects  measurements  to  produce  the 
topocentric  pointing  vectors  l first,  l mid, ,  and  liast.  Calculations  proceed  as  previously 
described. 
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Appendix  C.  Select  Satellite  Observations 


IN  the  course  of  three  sorties,  the  following  satellites  were  tracked.  Site  parameters  for 
all  observations  are  listed  in  Table  C. 


Table  C.l:  Site  Parameters  for  Select  Observations 


Geodetic  Latitude  </>: 
Geodetic  East  Longitude  A e 
Height  Above  Geoid  H 


39.6802°  (North) 
-83.8383°  (West) 
287.6  [m] 


Table  C.2:  Select  Satellite  Observations 


Satellite  Name  (NORAD  ID) 

Observation  Times  Azimuth  Elevation 

THOR  AGENA  D  R/B  (733) 

16- Jan-2008  10:59:00Z  220.76  60.26 

16-Jan-2008  10:59:38Z  211.08  46.35 

16- Jan-2008  11:00:19Z  206.61  34.78 

COSMOS  1812  (17295) 

2 8- Jan-2008  00:08:26Z  66.89  57.58 

2 8- Jan-2008  00:09:10Z  115.15  53.09 

2 8- Jan-2008  00:09:17Z  120.51  51.01 

COSMOS  1825  (17566) 

2 8- Jan-2008  00:24:37Z  321.66  47.38 

2 8- Jan-2008  00:25:21Z  340.87  34.07 

2 8- Jan-2008  00:26:05Z  349.80  24.02 

SL-14  R/B  (18215) 

03-Feb-2008  00:26:16Z  28.96  43.72 

03-Feb-2008  00:28:32Z  133.73  42.25 

03-Feb-2008  00:30:33Z  155.72  17.54 

SL-3  R/B  (19046) 

2 8- Jan-2008  00:11:11Z  162.70  68.06 

2 8- Jan-2008  00:12:42Z  348.73  50.64 

2 8- Jan-2008  00:13:59Z  347.90  25.78 

SL-16  R/B  (19120) 

16- Jan-2008  10:49:52Z  302.56  81.97 

16- Jan-2008  10:50:34Z  176.08  75.15 

16- Jan-2008  10:51:16Z  166.06  57.62 

COSMOS  1980  (19649) 

03-Feb-2008  00:37:13Z  172.87  45.03 

03-Feb-2008  00:39: 15Z  84.62  64.87 

03-Feb-2008  00:40:33Z  45.42  44.52 
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Table  C.3:  Select  Satellite  Observations,  Continued 


Satellite  Name  (NORAD  ID) 

Observation  Times  Azimuth  Elevation 

SL-8  R/B  (20433) 

16- Jan-2008  11:03:26Z  184.16  59.07 

16- Jan-2008  11:04:45Z  51.03  76.60 

16- Jan-2008  11:05:27Z  29.67  57.59 

ERS  1  (21574) 

2 8- Jan-2008  00:39:27Z  65.50  51.17 

2 8- Jan-2008  00:40:00Z  44.54  47.67 

2 8- Jan-2008  00:40:10Z  39.38  46.01 

UARS  (21701) 

03-Feb-2008  00:43:11Z  203.98  29.90 

03-Feb-2008  00:43:20Z  202.30  32.68 

03-Feb-2008  00:43:51Z  192.70  45.00 

COSMOS  2219  (22219) 

16-Jan-2008  10:38:19Z  91.94  55.35 

16- Jan-2008  10:38:55Z  70.96  49.89 

16- Jan-2008  10:39:03Z  67.28  48.29 

SL-14  R/B  (22287) 

16- Jan-2008  10:56:33Z  112.03  73.84 

16- Jan-2008  10:56:34Z  114.68  73.25 

16- Jan-2008  10:56:36Z  117.28  72.91 

SL-16  R/B  (23088) 

03-Feb-2008  00:11:45Z  70.01  29.44 

03-Feb-2008  00:12:44Z  86.91  25.91 

03-Feb-2008  00:13:27Z  96.60  22.49 

SL-16  R/B  (23705) 

16- Jan-2008  11:26:45Z  349.90  31.36 

16- Jan-2008  11:27:25Z  354.70  39.84 

16-Jan-2008  11:28:24Z  9.75  55.93 

COSMOS  2333  (24297) 

03-Feb-2008  00:46:23Z  136.29  43.73 

03-Feb-2008  00:46:31Z  2454499.53  132.46 

03-Feb-2008  00:47:06Z  114.45  47.87 

SL-8  R/B  (27535) 

2 8- Jan-2008  00:34:22Z  96.58  45.62 

2 8- Jan-2008  00:34:56Z  110.05  42.07 

2 8- Jan-2008  00:35:00Z  111.40  41.60 

SL-16  R/B  (28353) 

03-Feb-2008  01:06:54Z  171.70  47.70 

03-Feb-2008  01:07:03Z  169.18  50.16 

03-Feb-2008  01:07:37Z  153.93  60.02 

CZ-4B  R/B  (29093) 

03-Feb-2008  00:58:29Z  328.46  33.16 

03-Feb-2008  00:58:37Z  330.14  30.80 

03-Feb-2008  00:59:09Z  334.82  23.23 

SKYMED  1  (31598) 

16- Jan-2008  11:14:45Z  146.73  61.92 

16- Jan-2008  11:15:59Z  9.40  66.05 

16- Jan-2008  U:16:42Z  357.03  45.87 

SL-16  R/B  (31793) 

27-Jan-2008  23:55:57Z  40.12  53.91 

27- Jan-2008  23:57: 19Z  30.74  32.81 

27- Jan-2008  23:59: 16Z  26.94  15.96 
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